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ABSTRACT 



In this work we conclude the analysis of our CO line survey of Luminous In- 
frared Galaxies (LIRGs: Lir>10^^ L©) in the local Universe (Paper I), by focusing 
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on the influence of their average ISM properties on the total molecular gas mass 
estimates via the so-called Xco=M(H2)/Lco,i-o factor. One-phase radiative trans- 
fer models of the global CO Spectral Line Energy Distributions (SLEDs) yield an 
Xco distribution with: (Xco)~(0.6 ± 0.2) Mq (Kkms^-*^ pc^)^-*^ over a significant 
range of average gas densities, temperatures and dynamical states. The latter 
emerges as the most important parameter in determining Xco, with unbound 
states yielding low values and self-gravitating states the highest ones. Neverthe- 
less in many (U)LIRGs where available higher-J CO fines (J=3-2, 4-3, and/or 
J=6-5) or HCN fine data from the literature allow a separate assesment of the 
gas mass at high densities (>10'^cm~^) rather than a simple one-phase analysis 
we find that near- Galactic X co'^ (3-6j ]\/[q kTfi s ^ pc^ J ^ values heconne possi- 
ble. We further show that in the highly turbulent molecular gas in ULIRGs a 
high-density component will be common and can be massive enough for its high 
Xco to dominate the average value for the entire galaxy. Using solely low-J CO 
lines to constrain Xco in such environments (as it has been the practice up to 
now) may have thus resulted to systematic underestimates of molecular gas mass 
in ULIRGs as such lines are dominated by a warm, diffuse and unbound gas 
phase with low Xco but very little mass. Only well-sampled high-J CO SLEDs 
(J=3-2 and higher) and/or multi-J observations of heavy rotor molecules (e.g. 
HCN) can circumvent such a bias, and the latter type of observations may have 
actually provided early evidence of it in local ULIRGs. The only way that the 
global Xco of such systems could be significantly lower than Galactic is if the av- 
erage dynamic state of the dense gas is strongly gravitationally unbound. This is 
an unlikely possibility that must be nevertheless examined, with lines of rare iso- 
topologues of high gas density tracers (e.g. H^^CN, high-J ^^CO lines) being very 
valuable in yielding (along with the lines of the main isotopes) such constraints. 
For less IR luminous, disk-dominated systems, we find that the galaxy- averaged 
Xco deduced by one-phase models of global SLEDs can also underestimate the 
total molecular gas mass when much of it lies in a SF-quiescent phase extending 
beyond a central star-forming region. This is because such a phase (and its large 
Xco) remain inconspicious in global CO SLEDs. Finally detailed studies of a 
subsample of galaxies finds ULIRGs with large amounts (~10^ Mq) of very warm 
(>100K) anc? dense gas (>10^ cm~^) that could represent a serious challenge to 
photon-dominated regions as the main energy portals in the molecular ISM of 
such systems. 

Subject headings: galaxies: ISM — galaxies: starburst — galaxies: AGN — 
galaxies: IRAS — ISM: molecules — ISM: CO 
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1. Introduction 

Soon after the discovery of the luminous infrared galaxies (LIRGs), whose bolomet- 
ric luminosities are dominated by the infrared part of their Spectral Energy Distributions 
(SEDs) (LiR>lO^^L0) (e.g. Soifer et al. 1987), single dish and interferometric CO J=l-0, 
2-1 line observations were used to determine their total molecular gas mass and its distri- 
bution (Sanders et al. 1988a; Tinney et al. 1990; Wang et al. 1991; Sanders et al. 1991; 
Solomon et al. 1997; Downes & Solomon 1998; Bryant & Scoville 1996, 1999). These efforts 
were paralleled by several investigations of the so-called Xco=M(H2)/Lco(1^0) factor and its 
dependance on the average ISM conditions both theoretical (Dickman et al. 1988; Maloney 
& Black 1988; Wolfire et al. 1993; Sakamoto 1996; Bryant & Scoville 1996; Wall 2007) and 
observational (e.g. Israel 1988, 1993, 1997; Solomon et al. 1997; Downes & Solomon 1998; 
Yao et al. 2003). The average molecular gas conditions in LIRGs used in such studies have 
been typically constrained using CO(2-1)/(1-0) and CO/^^CO J=l-0, 2-1 line ratios (e.g. 
Braine & Combes 1992; Horellou et al 1995; Aalto et al. 1995; Papadopoulos & Seaquist 
1998). Higher- J transitions (J=3-2 and higher) were used only sporadically and mostly for 
star-forming galactic nuclei (e.g. Devereux et al. 1994; White et al. 1994; Giisten et al. 
1996; Nieten et al. 1999; Mauersberger et al. 1999; Dumke et al. 2001; Yao et al. 2003). 
This was a result of the larger difficulties such observations pose in terms of available submm 
receivers, their sensitivity, and the dry weather conditions needed (especially for z/>460 GHz, 
CO J=4-3). Receiver sensitivity limitations also hindered large multi-J line surveys of the 
much fainter lines from heavy rotor molecules such as HCN that probe higher density gas 
(>10^cm~^) except in nearby galactic nuclei (Jackson et al. 1995; Paghone et al. 1997) and 
a few luminous ULIRGs (e.g. Gracia-Carpio et al. 2008). 

Such limitations will soon be overcome after the ongoing commissioning of the Ata- 
cama Large MiUimeter Array (ALMA) is completed. Then routine multi-J observations of 
CO and heavy rotor molecules will yield unhindered view over the entire range of physical 
conditions in molecular clouds, from their quiescent and low-density phase (n(H2)~(10^- 
10^) cm~^, Tkin~(10-15) K) to the dense and warm gas intimately associated with star for- 
mation (n(H2)~(10^-10^) cm ^, Tkin~(30-150) K). The power of interferometric multi-J line 
imaging in revealing the mass distribution of dense warm SF gas in LIRGs has already been 
demonstrated by pioneering SMA observations (Sakamoto et al. 2008; Wilson et al. 2009; 
lono et al. 2007, 2009), while in the grand spiral M51 CO line ratio imaging at high res- 
olution revealed AGN-excited gas in its nucleus (lono et al. 2004). The influence of the 
high-excitation conditions found in SF regions gas on the Xco in galaxies may not necessarily 
be strong since dense and warm SF gas amounts to only ~(0.5-3)% of typical Giant Molec- 
ular Clouds (GMCs) mass. Even smaller fractions of the total molecular gas in spirals disks 
resides in their centers (~(0.1-1)%) where strong tidal fields, high cosmic ray energy densities 
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and/or AGN can drive a high molecular line excitation. Nevertheless this may no longer be 
true for the merger-driven starbursts in ULIRGs where a dense SF gas phase can contain 
the bulk of their total molecular gas mass (e.g. Solomon et al. 1992; Gao & Solomon 2004). 
Moreover, cases of AGN-driven mechanical and radiative feedback affecting the bulk of the 
molecular gas of the host galaxy and the corresponding CO SLEDs have now been identified 
(Papadopoulos et al. 2008; van der Werf et al. 2010). These systems along with ULIRGs, 
yield a nearby glimpse of ISM conditions that could be prevelaihng in the distant Universe. 

In the present work wc examine the influence of the average molecular gas conditions 
found in LIRGs (Papadopoulos ct al 2011, hereafter Paper I) on the Xco factor. We do so by 
using the largest combined database of LIRGs/CO transitions for which such a study has 
been conducted, while discussing also the limitations and potential biases of past theoretical 
and observational studies. We then outline methods that could be employed in the upcoming 
era of ALMA, and the special role the Herschel Space Observatory (HSO) can play, towards 
improved total molecular gas mass estimates, especially for ULIRGs (Lir>10^^ L©). Several 
such galaxies whose CO line ratios indicate extreme ISM conditions (see Paper I) are now 
studied individualy, their impact on the Xco values examined in detail. Throughout this 
paper we adopt a flat A-dominated cosmology with Ho=71kms~^ Mpc~^ and r2in=0.27. 



2. Moleculeir gas physical conditions and mass estimates in LIRGs 

The formal dependance of the Xco factor on the average density, temperature, and kine- 
matic state of large molecular cloud ensembles (where the statistical notion of Xco remains 
applicable) is explored in several papers (e.g. Dickman et al. 1986; Young & Scoville 1991; 
Bryant & Scoville 1996; Solomon et al. 1997; Papadopoulos & Seaquist 1999; Downes & 
Solomon 1998; Yao et al. 2003). CO and ^^CO lines can yield constraints on these ISM 
properties, and thus on the corresponding Xco, via radiative transfer models (e.g. Mao 
et al. 2000; Weiss et al. 2001). In this regard low- J CO SLEDs (up to J=3-2) with 
ncrit~(400cm~^-10^) cm""^ and Ej/Kb~ (5.5-33) K are adequate for determining the average 
state of the molecular gas and thus the appropriate Xco, provided that most of its mass is 
distributed in ordinary GMCs. The low-J CO SLED segment and its modeling can then in 
principle yield the mass normalization for the entire CO or of any other molecular line SLED 
(e.g. of HCN) typically emanating from much smaller mass fractions. For several ULIRGs 
in our sample whose CO SLEDs will be extended up to J=13-12 using the HSO this nor- 
malization is especially important as it allows determining the mass of the highly-excited 
molecular gas emitting the very high-J lines, and setting important constraints on the energy 
source responsible for its excitation (van der Werf et al. 2010; Rangwala et al. 2011). 
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2.1. Prior work and some methodological limitations 

There are currently only three major observational studies of Xco using CO lines for sub- 
stantial LIRG samples (Solomon et al. 1997; Downes & Solomon 1998; Yao et al. 2003). For 
the highly turbulent molecular gas in (U)LIRGs these typically find Xco=(l/10-l/3)Xco,Gai5 
with the low values attributed mostly to non self-gravitating gas distributions. For SF spiral 
disks a Xco,Gai~ (4-5) X0 remains applicable as most of their molecular gas is found in cool, 
low-density, self-gravitating GMCs as those in the Milky Way, pockmarked by SF "spots" of 
warm and dense gas with also nearly Galactic Xco (e.g. Young & Scoville 1991). In isolated 
spirals low Xco~(l/20-l/5)Xco,Gai values can be found only in their nuclei (2r<100-200pc) 
(Regan 2000), but involve only small fractions of the total molecular gas reservoirs. Thus for 
the metal-rich environments of LIRGs current observational work points to a bimodal Xco dis- 
tribution, with near-Galactic values for all isolated spiral disks, and (Xco)~(l/4-l/5) xXco,Gai 
in the merger-driven starbursts of ULIRGs. This view has even been widely adopted even 
for high-z star-forming LIRGs where only sparse sampling of their CO SLEDs exists (e.g 
Greve et al. 2005; Tacconi et al. 2006; Dannerbauer et al. 2009). It is worth noting that 
the claimed bimodality of the Schmidt-Kennicutt star formation relations (linking SFR and 
molecular gas supply) between disks and mergers (Daddi et al. 2010a; Genzel et al. 2010) is 
nearly equivalent to such an Xco bimodality (Narayanan et al. 2010 and references therein). 

Nevertheless much of the aforementioned observational work contains some serious 
methodological limitations borne out from the limited molecular line data available per 
object and specific assumptions made about the radiative transfer models. These are: 

• Setting Tkin=Tdust in the CO line radiative transfer models used to constrain Xco even 
though photoelectric, turbulent, or cosmic-ray heating can easily induce Tkin^Tjust- 

• Using a constant dV/dR (usually lkms~^ pc^^) in the Large Velocity Gradient (LVG) 
radiative transfer models that interpret CO line ratios even if: a) virial gas motions 
alone correspond to ~8-10 times higher values for the dense gas (n>10^ cm~^) that 
may be dominant in ULIRGs, and b) stellar mass concomitant with molecular gas 
and/or tidally- induced velocity fields can easily yield much larger dV/dR values. 

• The use of empirical relations and object-specific assumptions for obtaining Xco ex- 
pressions that can account for the potentially significant stellar mass concomitant with 
the molecular gas, an issue that arises especially in ULIRGs. 



1Xi=Mq (Kkms-ipc^)-! 
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• The molecular lines used in such studies (low-J CO lines) are insensitive to the prop- 
erties of the gas with densities n(H2)>10^ cm~^. 

The last two set the most serious methodological limitations on the very influential study 
of Xco in ULIRGs by Downes & Solomon 1998 (herafter DS98), whose results are widely 
used for presumably similar systems at high redshifts. In that study models of the CO 
J=l-0, 2-1 interferometric images of 5 ULIRGs were used to deduce Xco~(l/5-l/6)Xco,Gab 
attributing these low values to a continous molecular gas distribution encompassing large 
fractions of stellar mass and thus velocity fields determined by the total (gas) + (stellar) 
mass (see also Downes et al. 1993; Solomon et al. 1997). The use of empirical relations 
(e.g. LiR,/L^Q=2OOL0 (Kkms"-*^ pc^)""*^), and setting a rin_*=Mnew,*/Mbuige,*~0.5 mass ratio of 
newly formed stars (Mnew,*) to those in an old stellar bulge (Mbuige,*) in the pre-merger spirals 
make the computed Xco specific to the few ULIRGs studied and certainly not automatically 
applicable to other ULIRGs in the local or the distant Universe. Moreover "freezing" dV/dR 
to 1 kms^-*^ pc^-*^ renders such models incapable of constraining the average dynamical state 
of the gas and thus exploring its efi'ect on Xco in a straightforward way (see section 2.2). 

Furthermore, while the DS98 formalism can be generalized to arbitary r^ ,^ values, it 
remains impractical for the dusty gas-rich systems at high redshifts where stellar populations 
are heavily dust-enshrouded. Finally an Xco factor deduced solely from CO J=l-0, 2-1 line 
emission models may be inapplicable for a much denser gas phase (n(H2)>10'^ cm~^). This 
will yield a small error when such a phase represents only a few% of the total mass per 
typical CMC (as is the case in spiral disks), but to a potentially much larger error if the 
dense gas dominates the molecular gas mass budget. These difficulties are compounded 
by setting dV/dR=l kms~^ pc~^, which while justified by the sole CO(2-1)/(1-0) line ratio 
available in the DS98 study, may be partly responsible for the contradictory results obtained 
for the average gas density and CO line excitation. Indeed while HCN and high-J CO 
observations indicate much of the gas in some ULIRGs to be very dense (n(H2)>10'^ cm~"^), 
radiative transfer models of CO J=l-0 and J=2-l interferometric images of some of the 
same ULIRGs (e.g. Arp 220) become incompatible with these images already for n>10'^ cm~^ 
(DS98). Velocity gradients of dV/dR»l kms~^ pc~^ (expected for the highly turbulent gas 
disks of ULIRGs) and high temperatures for the dense gas can much reduce fine optical 
depths and the emergent CO line emission possibly aleviating these discrepancies. 

In summary all current observational studies of Xco in LIRGs, leave its range and de- 
pendance on the average density, temperature and kinematic state of the molecular gas 
still largely unexplored. The small LIRG/CO line datasets used in such studies, and the 
aforementioned limitations of the analysis methods, limits their applicability to other such 
systems in the local or the distant Universe. 
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2.2. The average state of the moleculeir gcis and Xco 

Our dataset of total CO lines luminosities for a substantial sample of LIRGs provides a 
good opportunity for exploring the Xco dependance on the average ISM state in vigorously 
star-forming systems. The average dV/dR for the gas is now a free parameter to be con- 
strained by LVG radiative transfer models along with the density and temperature. This 
provides an extinction-free and straightforward method for determining the effects of stellar 
mass concomitant with molecular gas and/or tidal disruption of molecular clouds on Xco. 

To place our discussion regarding molecular gas mass estimates in LIRGs in an au- 
tonomous context we reproduce the Xco factor (see Appendix) as 



^ M(H2) ^ 3.25 VK(H^ 1 f \ 
°° Lco(i_o) Tb,i-o ^^Kkms-ipcV' 

where q;=0.55-2.4 depending on the assumed cloud density profile (see Bryant & Scoville 
1996, Eqs A 11, A 17), and n(H2), Tb,i-o are the average gas density and CO J=l-0 brightness 
temperature of the molecular cloud ensemble. The parameter Kyir is given by 



(dV/dR)^^^ VlO'cm 

(with Aco = [G0/H2]/(dV/dR being one of the outputs of a typical LVG model), and deter- 
mines the average dynamic state of the molecular gas, with Kvir~l-3 for the (mostly) self- 
gravitating GMCs in the Galactic disk (and Kvir^l corresponding to dynamically unattain- 
able gas motions). For a typical value of a=1.5 Equation 1 becomes 



Xco = 2.65^-^K7i ^ , (3) 

Tb,i-o VKkms ipcV 

which we use in the present work. For ordinary GMCs with n(H2)~500cm~^, Tb,i_o~10K, 
and Kvir~l-2 the latter yields Xco~(3-6)Xi (where Xi=M0(Kkms-^ pc^)-!). 

The optically thin limit for CO J=l-0 line emission is often adopted to provide a lower 
limit on the total molecular gas mass (e.g. Bryant & Scoville 1996; Solomon et al. 1997). 
The availabihty of several CO lines allows a more robust such hmit without assuming LTE: 
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Xll^'"^) = 0.078 



. 1 55/T ,n 5 7 2J + 3 
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K km s ^ pc^ 



(4) 



(see Appendix), which makes obvious that, unless the CO hne ratios become significantly 
larger than unity (which they can in optically thin emission), most of the contributions come 
from the low J levels. The LTE expression is: 



which gives a higher (and less reliable) lower limit since the LTE partition function can be 
significantly larger than the non-LTE one. The optically thin approximation can be used also 
for the realistic case of significant CO line optical depths as long as large velocity gradients 
keep the photon escape probability local throughout the bulk of the molecular gas mass. 
Then Equation 4 is modified to 



K ' 



Ma 







Kkms ^ pc^y ' 
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(6) 



Kkms ^pc^ 



where /3j+ij=[l — exp(— rj_|_ij)] /rj+u is the photon escape probability (see Appendix). The 
last equation provides a more robust lower limit than Equation 4, as it accounts for finite line 
optical depths (computed from radiative transfer models). If enough line ratios arc available, 
the estimate from Equation 6 will approach that of Equation 3 as long as one average gas 
phase dominates the emergent CO line emission. 

In Figure 1 we show the Xco and Xco'^"^'' (LTE) distributions, computed from Equations 3 
and 5 and the results of one radiative transfer models for the CO SLEDs and ^"^CO lines for 
all the LIRGs in our sample (see Paper 1). Most Xco values are <L5Xi, with (Xco)~0.60Xi 
obtained for the main distribution (excluding the few outliers beyond 3Xi). The Xco'^™'' (LTE) 
distribution gives similar results, indicating that most of the states compatible with the CO 
SLEDs have optically thin or moderately optically thick CO J=l-0 hne (rio~l-3). This has 
been noted in the past in the context of a two-phase ISM model and attributed to the high 
temperatures of a turbulent "envelope" phase in molecular clouds (Aalto et al. 1995). 

The Xco range obtained using our one-phase models of global CO SLEDs is similar 
to that reported by DS98 and Yao et al. 2003. Moreover in Figure 2 the distribution of 
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y^n(H^/Tb,i_o demonstrates that unlike often stated, the effects of density and temperature 
do not always cancel out but arc responsible for much of the Xco variations over different 
ISM environments (compare Figure 2, with the Xco distribution in Figure 1). Nevertheless 
the average dynamical state of the gas (i.e. Kvir) remains the most important factor affecting 
the average Xcq. This is shown in Figure 3 where the Xco distributions for virial/near-virial 
and unbound average gas dynamical states are shown. In the merger-driven starbursts of 
ULIRGs highly non-virial gas velocities can be caused by significant amounts of stellar mass 
concomitant with a continous molecular gas distribution as noted by DS98 (see also Solomon 
et al. 1997), albeit without using radiative transfer models to constrain Kvir- Strong tidal 
fields acting on GMC cloud envelopes and/or a diffuse intercloud molecular gas distribution 
whose velcity field traces the combined potential of the surviving dense GMC cores as well as 
stars will both raise the average Kyir, and further reduce the corresponding Xco (Equation 3). 

In Figure 4 the distribution of Xco^ computed from Equation 6 is shown which, unlike 
that deduced from Equation 3, uses the CO line ratios explicitely. This makes it more 
suitable when fully-sampled CO SLEDs from J=l-0 up to high-J (e.g. J=6-5, 7-6) are 
available, and when the turbulent models of the hierarchical density structures of molecular 
clouds (e.g. Ossenkopf 2002) become refined enough to provide global CO and ^"^CO SLEDs 
as well as the corresponding values per density "sub-phase". For well-sampled ^^CO 

SLEDs the X^^^ expression (using the appropriate [^^C0/H2] abundance) is even more useful 
as ^j+i,j(^^CO)~l. In our current study using the output values from our one-phase 

models in Equation 6 yields an Xco^ distribution similar to that of Xco from Equation 3. 

Interestingly the Xco distributions, while concentrated within ~(0.3-L5) Xi, do extend 
out to Galactic values ~(2.5-6) Xi. These are found for disk-dominated LIRGs (e.g. NGC 157 
and Mrk 1048), (U) LIRGs with unexpectedly "cold" CO ratios (e.g. IRAS 05189-2524) where 
cool (~15 K) gas with low/modest densities (~(10^-10^) cm~^) dominates their average ISM 
state. Surprisingly, Galactic X^o values can he found also in some ULIRGs, whose "hot" 
CO ratios (e.g. IRAS 08572-1-3915) indicate a dominant warm (~(100-150) K) and dense 
(~(3x 10^-10^) cm~"^) phase, or other evidence (e.g. multi-J HCN lines) suggest a massive 
dense, gas phase (e.g. Mrk 231). Thus large Xco values are possible for both high and low 
excitation gas phases and thus cannot be used as indicators of Galaxy-type ISM conditions 
in distant star forming disks (e.g. Daddi et al. 2010b). This also suggests that in many local 
ULIRGs and similar systems at high redshifts the widely adopted Xco^ (0.8-1) Xi values 
understimate their molecular gas mass. 
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2.3. Eddington-limited stcir formation and a minimum moleculcir gas mass 

Recent studies of SF feedback suggest a maximum eg ,,=Li^VM*(H2)~500 (L0/M0) for 
the dense and warm gas M*(H2) near SF sites in galaxies as a result of strong radiation 
pressure from the nascent O, B star clusters onto the concomitant dust of the accreted gas 
fueling these sites (Scoville 2004; Thompson et al. 2005; Thompson 2009). Thus, provided 
that average dust properties (e.g. its effective radiative absorption coefficient per unit mass) 
remain similar in metal-rich star-forming systems such as LIRGs, a near-constant * is 
expected among them. A value of eg^*~500 (Lq/M©) is actually measured in individual SF 
sites of spiral disks such as M51 and entire starbursts such as Arp220 (Scoville 2003), while 
~(440 ± 100) (Lq/Mq) is obtained for CS-bright star-forming cores in the Galaxy (Shirley 
et al. 2003). Further evidence for a near-constant Cg .^ is the tight linear Lir-HCN(I-O) 
correlation found for individual GMCs up to entire ULIRGs (Wu et al. 2005) with HCN 
J=l-0 used as a dense gas mass tracer. However the intermitency expected for galaxy-sized 
molecular gas reservoirs (i.e. at any given epoch of a galaxy's evolution some dense gas 
regions will be forming stars while others will not) will lower the global eg^* to ~ 1/3- 1/2 of 
the Eddington value (Andrews & Thompson 2011). 

Similar eg,* values can be obtained without explicit use of the Eddington limit (and the 
detailed dust properties it entails), but from the typical L*/Mnew,* in young starbursts where 
Mnew,* is the mass of the new stars and L* their bolometric luminosity ( for the deeply 
dust-enshrouded SF sites). For esF,c=Mnew,*/[Mnew,* + M*(H2)] as the SF efficiency of the 
dense gas regions where the new stars form it is 



For esF,c~0.3-0.5 typical for dense SF regions, and Lj*^/M„ew,*=(300-400) (Mq/Lq) (Downes 
& Solomon 1998 and references therein). Equation 7 yields eg,^,~( 130-400) (M0/L0). Here 
we choose eg,*=250 (Lq/M©), close to the average values from Equation 7 and the black body 
limit deduced for compact CO line emission concomitant with an optically thick (rioo/im>l) 
dust emission (Solomon et al. 1997). Eddington-limited star formation in LIRGs imphes a 
minimum molecular gas mass MsF—Lm/^g,* fueling their observed star formation rates. In 
Figure 5 we show the Msf distribution which, for the ULIRGs in our sample, reaches up to 
~5xlO^M0 i.e. surpassing the total molecular gas reservoirs of typical spirals. 




(7) 
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2.4. Computing Xco in a two-component approximation 

In some cases where one-phase LVG models of global CO line ratios of LIRGs do not 
converge to a well-defined range of physical conditions, a simple two-phase model can be 
used to examine the underlying ISM conditions and the corresponding Xco factors. The 
total CO J=l-0 luminosity can then be expressed as 



■'-'CO, 1-0 — ■'-'00,1-0 + ■'-'00,1-0 ~ l^g,*^co ) ^IR ■'■^00,1-0' \°) 

where (h) and (1) indicate the high and the low excitation phases, Xco^ is the CO-H2 con- 
version factor for the former, and eg,*=25O(L0/M0). This assumes that the high excitation 
phase fuels Eddington-limited star formation. For any other transition 



-■CO.J+l-J 



L 



(*) 

IR 



tgj*^CO 



Jh) ,t(1)' 
^J+IJ ^ ■'-'00,J+1-J' 



(9) 



The molecular gas near HII regions and interfaces between molecular clouds and 0,B 
stellar associations is an obvious choice for obtaining template high-excitation CO SLEDs. 
The beam- matched CO line survey of the Orion A and B GMCs found r2i=1.2-1.3 and 
R2i=10 (the CO/^^CO J=2-l ratio) in such "hot" spots (Sakamoto et al. 1994). These 
ratios along with Tkin>100K (expected from the high CO line brightness temperatures 
in such SF spots) set as constraints to a one-phase LVG model yield: Tkin=(125-150) K, 



n(H2)=3xlO^ cm ^, Kvir=7, and corresponding line ratios: r^'j^(3)=1.35, r^2''=l-33, r^3''=1.30. 



.(h). 



.(h) _ 



,(h). 
'■54 



,(h)_ 
'65 



=1.25, and rlp'' = 1.22. For this LVG solution set we compute (Equation 3) Xco^ 
~2.2Xi, and from Equations 8 and 9 the CO line ratios of the low-excitation (l)-phase. 
These can then be used as constraints on 1-phase LVG models to obtain the average ISM 
conditions and Xco • For such a 2-phase decomposition the effective xio ^'^^^ factor would be 



y(h) ^(l-h)y(l) 
Y(2-ph) _ ^co fJco ^co 



(10) 



with pco '^^=Lco i-o/Lco^-o- Using the Eddington- limit normalization for Lco\_o yields 



(h)' 



(h)' 



T (*) / \ 

x^^^) = x« + 1 - ^ I (11) 

%*^oo,i-o \ X^o / 

When lines such as higher-J CO transitions (J=4-3, 6-5) and/or heavy rotor molecular lines 
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(e.g. HCN) are available, they are used to set constraints on the (h)-phase, pco , and 
Xco"'''^'' without assuming the (h)-phase CO SLED and X^o^ deduced from the Orion CMC 
star-forming "spots" . 



3. The Xco factor in LIRGs 

The average ISM conditions deduced using the CO lines available per LIRG (see Paper 
I) encompasses the [n(H2), Tkin, Kyir] parameter space where most of the mass is expected to 
reside in ordinary GMCs. Thus the narrow Xco distribution around (Xco) ~0. 6X1 (Figs 1, 4), 
does not seem as the result of an excitation bias induced by the particular molecular lines 
used. This is not to say that other high- J CO or e.g. HCN line luminosities are expected to 
be compatible with the average ISM states deduced from our current CO line dataset but 
rather than that such lines with their higher critical densities will probe much smaller gas 
mass fractions per typical CMC to be of any consequence when it comes to the global Xco- 

It is tempting to consider the aforementioned picture as complete. It is certainly com- 
patible with models of supersonic turbulence in ordinary GMCs where most molecular gas 
mass is found at densities <10^cm~^, its conditions thus "accessible" to the CO, ^'^CO lines 
used to constrain them (and Xco) in our study. In the merger-driven starbursts of ULIRGs 
however GMCs may be far from ordinary. Stripping of their outer envelopes by strong tidal 
fields and/or the large pressures from a hot ionized gas phase resulting from HI cloud col- 
lisions, both expected in merger environments (e.g. Solomon et al. 1997), can dramatically 
alter the M(H2)-n(H2) distribution towards most mass being at n(H2)>10^ cm~^. This ac- 
tually occurs in the Galactic Center (e.g. Giisten & Philipp 1994), and in ULIRGs it can 
involve their entire molecular gas reservoirs. In such cases low-J CO, ^^CO lines are no longer 
sensitive to the average ISM conditions and may thus yield inaccurate "corrections" to the 
global Xco factor with respect to its Galactically calibrated value. 



3.1. The Xco in ULIRGs: theoretical expectations for highly turbulent gas 

The high pressures in the highly turbulent gas disks in ULIRGs (DS98) will: a) result 
to much larger average gas densities at all scales, and b) "relocate" large mass fractions of 
their GMCs to densities n(H)>10'^ cm""^. The latter can be easily shown from the probability 
distribution function (pdf) of the density in supersonically turbulent clouds. This is well- 
approximated by a log-normal distribution with a dispersion: 



3M2 



) 



nl/2 



In 1 + 



4 



(12) 



where M=(Tv/cs is the average 1-dim Mach number (Padoan & Nordlund 2002). The mass 
fraction contained in cloud structures with overdensities x>Xo (x=n/(n)) is then given by 



The high velocity dispersions (7v~ (30-140) km s~^ measured in the molecular disks of 
ULIRGs (DS98; see also Swinbank et al. 2011 for a recently discovered such disk at z~2.3) 
versus those in spiral disks (~(5-10) kms~^) correspond to M(ULIRGs)~(3-30) xM(spirals). 
For M(spiral)=10 then (Tp(spirals)~2, while (Tp(ULIRGs)~2.55-3.33, which dramatically ex- 
tends the gas density pdf expected in ULIRGs towards high values where much of the gas 
mass will now lie. Indeed for the GMCs in spiral disks only ~3% of the mass will be found 
at overdensities of x>Xo=500 (=5xl0^cm^^ for typical GMC with (n)=100cm~^) while 
for ULIRGs this is ~45% (Figure 6). This difference can be even larger since the average 
molecular gas density in disks is lower ((n)~(100-500) cm~^) than in ULIRGs where it can 
reach up to at least ~10^cm~^ (Greve et al. 2009). Thus gas at n>5xl0^cm~^ corresponds 
to an overdensity of Xo=100-500 (containing <10% of gas mass) for GMCs in spirals, and 
only Xo=5 (and 90% of the mass) for GMCs in ULIRGs. The typical CO SLED available 
for LIRGs in our sample (J=l-0, 2-1, 3-2) would then be insensitive to the state of the 
molecular gas in ULIRGs and thus unable to constrain the corresponding Xco factor, while 
it would remain adequate for such a task for the molecular gas in isolated spirals. 

Interestingly the few cases of near-Galactic Xco in ULIRGs were found when either 
higher-J CO lines were available or when a dense and warm gas phase was massive enough to 
be obvious even in low- J CO SLEDs (Paper I, and section 4). The Xco for such dense gas can 
be high, approaching and even surpassing Galactic values. Indicatively, for n=5xl0^cm~^ 
and warm gas with a thermalized CO J=l-0 line Tb,i_o~Tkin~(100-150) K, it is: Xco~(4- 
6)K;^iJXi (Eq. 3), which for self-gravitating gas (Kvir=l) corresponds to Galactic values. 

Moreover, for the highly turbulent, high-pressure ISM environments in ULIRGs the 
notion of molecular gas reservoirs reducible to ensembles of dicrctc GMCs may no longer 
apply, with much of the gas in the pre-merger GMCs redistributed in continous disks of 
~(100-300)pc diameter (DS98; Sakamoto et al. 2008). A warm, low-density, molecular 
gas phase, with large Kyir (>10) can then be generated, alongside the denser one, by the 
tidal disruption of GMC "envelopes" in merger environments, the high pressures turning 
intercloud CNM HI into H2, and the effects of supersonic turbulence driven at the largest 
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scales (e.g. Ossenkopf 2002). Such a diffuse phase, even if containing httle mass, can easily 
dominate the emergent global CO J=l-0, 2-1 line emission of ULIRGs, yielding "cold" 
CO(2-1)/(1-0) ratios and low X^o factors (a result of its low densities, high temperatures, 
and large K^[^ values). These low Xco values, while appropriate for the diffuse gas mass, may 
not be so for the bulk of the molecular gas that now resides at much higher densities. 

Thus in the very turbulent molecular gas reservoirs of ULIRGs a combination of: a) 
most of their mass residing at high densities (with potentialy large X^o factors), and b) the 
existence of a diffuse, warm, and unbound phase with little mass (and low Xco) dominating 
low- J CO SLEDs (and the average Xco determined from them), can easily result to a sys- 
tematic underestimate of their total molecular gas mass. All current observational studies 
likely suffer from this bias, expected to be most prominent in ULIRGs with high molecular 
gas surface densities. This includes the present study as mostly low-J CO SLEDs (up to 
J=3-2) are typically available per LIRG. The only exceptions are the few ULIRGs for which, 
a highly excited CO J=3-2 line, CO J=4-3, 6-5 lines, and/or HCN lines from the literature 
allowed the revealing of their massive dense gas reservoirs (see Section 4). Studies of Xco in 
similar systems at high redshifts (e.g. Tacconi et al. 2008) will be similarly affected. Only 
molecular SLEDs with high critical densities can overcome such an Xco bias in ULIRGs and 
see how many actually lie on the high end of the Xco distributions shown in Figures 1 and 4. 



3.1.1. Numerical simulations and the Xco factor: important caveats, and ways forward 

Recently GMC-sized numerical simulations that include H2 and CO formation/destruction 
and radiative transfer in MHD turbulent models where used to explore Xco and its depen- 
dance on average ISM conditions (Shctty ct al. 2011a,b). These simulations also include 
chemistry and thermodynamical effects but not SF-driven heating (via photons or CRs) or 
turbulent heating which could dominate in ULIRGs. Currently such studies explore Galactic- 
type GMCs whose Xco values they found similar to the one observed. Extending such numer- 
ical models to ISM conditions expected in ULIRGs is not straightforward as GMC boundary 
conditions such as ambient radiation fields, surface pressure, external gravitational field and 
its tidal terms (important for galactic centers of spirals and ULIRGs) are set constant in 
GMC-sized simulations or, in the case of tidal fields, are omitted altogether. Nevertheless in 
any realistic setting of GMCs in galaxies these boundary conditions do change, sometimes 
faster than internal cloud evolutionary timescales, and especially so in ULIRGs. 

Galaxy-sized numerical studies of Xco in disks and mergers recently presented by Narayanan 
et al. 2011 address the aforementioned important point of coupled GMC-galaxy evolution 
by tracking the evolution of GMC boundary conditions, but do so by necessarily adopting 
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some subresolution methods to follow in-cloud physics and keep the problem computation- 
ally tractable. They recover Galactic X^o values for disks, but ~5-10 times lower ones in 
mergers, with a large spread that makes a single ULIRG-appropriate X^o factor impractical. 
Nevertheless, for ULlRGs in particular, these simulations do not have (and could not have) 
the resolution necessary to track gas with n>10^ cm~^ (where most of their gas mass resides). 
Resolving the all-important kinematic state of such a dense gas component (self-gravitating 
or not?) within their compact gas reservoirs (~(100-300) pc) is currently impossible with 
galaxy-sized simulations, as is the explicit tracking of the large turbulence-regulated range 
of densities expected in ULIRGs. There the average densities can surpass ~10^cm~^, while 
the high levels of turbulence will maintain >60% of the mass at n>10^cm~^ (Figure 6). 
With the Tkin values computed in such simulations likely to be lower limits (as CR and 
mechanical SNR-driven heating are omitted), gas densities and the corresponding kinematic 
states remain the only parameters affecting the average Xco in metal-rich environments that 
can push its values either way. Indeed while additional heating can only raise Tkm and 
thus the CO J=l-0 brightness temperature (reducing Xco), high densities (>10'^cm^^) and 
self-gravitating states (Kvir~l) will act to raise Xco. 

GMC-sized simulations of the much more turbulent and denser molecular gas in ULIRGs 
are necessary for computing the corresponding Xco but these must now use new initial con- 
ditions with: a) much higher volumed-averaged gas densities (current ones use ~150cm~^), 
b) higher velocity dispersions (current ones use (Tv~2.4kms~^), and c) higher background 
CR energy densities (and thus non-neghgible CR gas heating). Such simulations and their 
emergent line intensities can be compared to much richer molecular SLEDs available now 
than in the past. These include the typical CO lines from J=l-0 up to J=4-3 observed from 
the ground, rising up to J=13-12 using the Hcrschcl Space Observatory (van der Werf et 
al. 2010), further complemented by multi-J observations of heavy rotor molecules like HCN, 
HCO^, CS (e.g. Grcvc ct al. 2009) as sensitivities vastly improve in the era of ALMA. 
Thus, rather than simply trust the Xco yielded by GMC-sized simulations as is currently 
done, these simulations will yield Xco values as a by-product of the much more difficult task 
of reproducing relative strengths of molecular line luminosities with critical densities ranging 
from ~100cm~^ up to ~10^cm~^. These encompass the entire range of hierarchical struc- 
tures in turbulent GMCs, from a highly turbulent low-density possibly non self-gravitating 
"CMC-envelope" phase at large scales up to self-gravitating compact dense gas cores with 
dissipated supersonic turbulence (e.g. Ossenkopf 2002). GMC-sized numerical simulations 
exploring a dense grid of ULIRG-type CMC boundary conditions and determining the corre- 
sponding Xco grid can then be used to inform galaxy-sized numerical simulations where such 
boundary conditions are tracked, while resolving the hierarchical structures of individual 
GMCs (and thus determining f(xo) from the simulations) remains out of reach. 
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3.2. High X CO in ULIRGs: past observational evidence ignored? 

As discussed in the previous sections, by relying on low-J CO SLEDs, most observa- 
tional studies of Xco may have imparted a bias towards low values in ULIRGs. Nevertheless 
past observations of heavy rotor molecules such as HCN J=l-0 have given early hints of a 
large scale M(H2)-n(H2) redistribution in these merger systems, which placed large fractions 
(>50%) of their total molecular gas mass at much higher densities (>10^cm~^) than in iso- 
lated gas-rich disk systems (e.g. Solomon et al. 1992; Gao & Solomon 2004). However a 
single transition of a heavy rotor molecule cannot set constraints on the gas temperature and 
density (let alone on its average dynamical state) useful enough to constrain the correspond- 
ing Xmoi factor well. Thus these early estimates of the dense molecular gas mass in ULIRGs 
are necessarily very uncertain as they relied on an Xjjcn factor derived for n~nc;i.it [HCN(l-O)], 
self-gravitating gas, and HCN J=l-0 brightness temperatures typical of the IR color tem- 
peratures of SF galaxies (e.g. Gao & Solomon 2004). Indicatively for n(H2)=5xlO^ cm~^, 
Tb,i_o(HCN)=40 K and K^ir^l Equation 3 apphed for HCN J=l-0 yields Xhcn=15 Xj, which 
can be higher still for higher gas densities, as indicated by multi-J HCN, CS and HCO+ hne 
observations of Arp220 and NGC6240 and where Xhcn~(20-35) Xj (Greve et al. 2009). 

For such high Xhcn values and L^cn,i-o/Lco,i-o=1/4-1/6 observed in ULIRGs (e.g. 
Solomon et al. 1992; Gao & Solomon 2004; Gracia-Burillo et al. 2012) it is easy to show 
that X^o"^*"^ (Equation 10) will yield near-Galactic or higher values in such systems. Indeed, 
assuming that the only contribution to the (h)-phase in Equation 10 comes from the HCN- 
bright gas (i.e. the HCN emission from the low-excitation (1) phase is neglible) we can 
determine pco '''' from 



J I (obs) 
(1-h) _ ^co,l-0 _ _ ^CO/HCN _ , . 

Pco (h)' (h) ' ^ ' 

^co,l-0 ^CO/HCN 

where rco/HCN~-'^co,i-o/LHCN,i-o ™- ULIRGs(=4-6), and rj^Q^y^^p^ that of the (h)-phase. The 
latter can be found using a radiative transfer model for densities typical of HCN-bright gas, 
an assumed gas temperature, and setting K^'j'j,^ = l. For n=10^cm~^, and Tkin=100K using 
an LVG model for HCN and CO lines we compute Tb,HCNi-o=41 K and t[|'^qj^_q=96 K. Thus 
rco/HCN=2-34, while Xhcn~20Xi and X^o^~8.6Xi, with pio"''^~0.71-1.56. Setting xS=0.5Xi 
(Figs. 1, 4) for the low-excitation (CO-bright but HCN-dark) phase, we obtain Xco~'^'^'*~(4- 
5) Xi, i.e. Galactic values. These are now due to much higher fractions of dense HCN-bright 
gas than in the Galaxy, with large Xhcn (and Xco^). If such Xco values are the norm in 
the high-pressured turbulent gas disks of ULIRGs, their dynamical masses would he domi- 
nated by the molecular gas, quite unlike isolated SF spirals. Indeed for all the ULIRGs in the 
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DS98 study whose dynamical mass is reasonably well-constrained (with both CO 1-0 and 2-1 

imaging) Galactic Xco values would yield Mgas/Mdyn~0.6-1, and compatible with the obvious 
limitation of Mgas/Mdyn<l within the observational uncertainties of Mdyn estimates. Regard- 
ing the latter it is is worth pointing out that low S/N for extended CO line emission regions 
in current interferometric maps as well as dynamically unrelaxed, non-circular, molecular 
gas motions (expected in strong mergers) typically cause underestimates rather than over- 
estimates of the actual M^yn- Surprisingly such a bias was deduced even for gas-rich but 
otherwise isolated disks with numerical simulations showing an observationally-determined 
Mdyn to be typically be ~30% less than the actual value (Daddi et al. 2010b). Needless to 
say this state of affairs becomes worse still of CO-imaged systems at high redshifts. 

Multi-J observations of heavy rotor molecules and/or high- J CO hues in ULIRGs can 
produce dense gas mass estimates via a (radiative transfer model)-constrained Xmoi factor. 
However the few such studies available (Papadopoulos 2007; Papadopoulos et al. 2007; Krips 
et al. 2008; Gracia-Carpio et al. 2008; Greve et al. 2009) are a testimony of how hard it 
is to overcome the degeneracies inherent in radiative transfer models of optically thick line 
emission from heavy rotor molecules. The n-Tkin solution degeneracies of the Large Velocity 
Gradient (LVG) radiative transfer models typically used in such studies translate to a wide 
range of Xhcn factors, ranging from ~10Xi and reaching up to ~(50-90)Xi (Krips et al. 
2008; Gracia-Carpio et al. 2008). Telhngly even the two studies with the largest number 
of dense gas tracer hues available per ULIRG, that of Mrk231 (Papadopoulos et al. 2007) 
and of Arp220, NGC6240 (Greve et al. 2009) cannot overcome such degeneracies, and yield 
LVG solutions that produce Xhcn~(10-20) Xi and ~(17-37) Xi respectively. Nevertheless 
most Xhcn values deduced in such studies are high enough to yield dominant amounts of 
dense gas mass in ULIRGs, with a near-Galactic or even larger than Galactic Xco 

This apparent contradiction of the low Xco factor advocated for ULIRGs by past low-J 
CO line studies (e.g. Solomon et al. 1997; DS98; Yao et al. 2003), and the larger effective 
Xco "^^^ implied by heavy rotor molecular line emission (mostly HCN) in such galaxies has not 
been much noted in the literature. In the rare cases where this discrepancy was discussed, 
the argument went for a downward revision of Xhcn by ~l/5 rather than a boost of the 
widely adopted ULIRG- value of Xco~lXi (e.g. Papadopoulos et al. 2007). It must be noted 
however that the possibihty of HCN-based molecular gas mass estimates raising the effective 
Xco in ULIRGs is mentioned in the seminal interferometric work by DS98. 
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3.2.1. What could lower Xco in ULIRGs 

Lower Xhcn values than those reported in existing studies are possible if higher gas 
temperatures and/or unbound gas motions (Kvir>l) prevail for the bulk of the dense gas in 
ULIRGs. In the most detailed such study Greve et al. (2009) do find warm and non-virial 
LVG solutions for the HCN-bright gas in Arp220 and NGC6240 with n(H2)=3xl05cm-3, 
Tkin=(45-120) K (Arp220) and Tki„=(60-120) K (NGC6240) while Kvir~2, 5 for Arp220 
and NGC6240 respectively. Choosing Tkin=100K as indicative for both systems (for a 
given density and Kvu- high Tkm's yield small Xmoi's), we find Tb,HCNi-o(Arp 220)~60 K, 
TKoi-o(Arp 220)~96 K and Tb,HCNi-o(NGC 6240)~40 K, tJ^;^^i_o(NGC 6240)~86 K (for their 
corresponding Kvir values). Then from Equation 3: Xhcn~12 Xi, Xco^~7.6 Xi for Arp 220 and 
Xhcn~7.3Xi, x£^~3.4Xi for NGC 6240. Thus using r[.°o/^HCN=5-9(Arp 220), 12.5(NGC 6240), 

and the LVG-computed: r2'(^/jjcj^~1.6(Arp 220), 2.15(NGC 6240), we find pilr''^~2.68(Arp 220), 
4.81((NGC 6240). Then from Equation 10 and the aforementioned numbers (and setting 
xS=0.5Xi) it is: X^o"P^^(Arp220)~2.4Xi and X^o"P''^(NGC6240)~1Xi. 

In the aforementioned example the warm "end" of a degenerate range of LVG solutions 
to a multi-J heavy rotor molecular line dataset yields an Xco~(l/5)Xco,Gai for NGC 6240 
while reducing that of Arp 220 to ~(l/2)Xco,Gai- Nevertheless such sets of solutions, with 
Kvir>l and Tkin significantly higher than T^ust are not typically considered optimal in most 
studies. Indeed Kvir~l is often used as a constraint on the LVG solution range of molecular 
lines tracing dense gas (e.g. HCN, CS) while solutions with Tkin~Tdust are considered as 
more appropriate for such a gas component since gas-dust thermal coupling is expected to be 
strong (e.g. see the dense gas model in Mrk231 by Papadopoulos et al. 2007). On the other 
hand recent theoretical and observational work has shown that in ULIRGs large temperatures 
with Tkin>Tdust are possible even with most of the gas mass at high densities (n>10^ cm^^) 
because of dominant CR and/or turbulent heating, both powered by the large SNR number 
densities in these galaxies (Papadopoulos 2010; Paper I; Rangwala et al. 2011). Nevertheless, 
while high Tkin certainly contributes in lowering Xco'* of the dense gas phase, it is the high 
Kyir values that are mostly responsible for this, a result noted also by previous studies 
(e.g. DS98). Indeed, setting Ky^^—l for the aforementioned LVG solutions for the dense gas 
phase in Arp 220 and NGC 6240 brings their X^o"^*"^ factors back up to ~4.5 Xi(Arp 220) and 
~3.3Xi(NGC6240). 

Clearly good constraints on the dynamical state of the dense gas in ULIRGs is of ut- 
most importance since while Kyir^l may be a safe assumption for the SF- fueling dense gas 
component of individual GMCs in typical spirals, it may not be so for the extreme ISM 
environments of ULIRGs. Strong tidal fields and/or "bottom" -stirred ISM by strong SF 
feedback (via multiple SNR shocks and the radiative feedback implied by Eddington-limit- 
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regulated SFRs) could induce Kvir>l even for the dense gas in such systems (see Paper I), 
lowering Xhcn and hence the corresponding Xco~^'^^ factor. 

3.3. Cold SF-quiescent molecular gas and large Xco in LIRGs 

Gas- rich disks containing a starburst in their central ~(l-2) kpc, but also large amounts 
of SF-quiescent ISM at larger galactocentric radii are often found in LIRGs with Lir<1O^^L0, 
and such configurations are even expected for weakly perturbed systems (e.g. Maiolino et al. 
1997 and references therein). Large amounts of cold (~(10-15) K) molecular gas beyond a 
warm component (~(35-40) K) confined within a star-forming central ~1 kpc has been shown 
as a general feature of LIRGs using CO J=2-l, 1-0 lines (Papadopoulos & Seaquist 1998). 
Sensitive dust submm continuum and CO J=l-0 imaging in individual LIRGs found cold 
dust and concomitant molecular gas with low CO J=l-0 brightness extending out to radii 
of at least ~3kpc (Papadopoulos & Seaquist 1999; Papadopoulos & Allen 2000). Submm 
continuum and HI imaging studies of nearby spirals also suggest cold molecular gas as a 
general feature of their disks at large galactocentric distances (Thomas et al. 2004). 

Extended, SF-quiescent, molecular gas reservoirs in the disks of LIRGs can be unde- 
tected because of low CO J=l-0 hue brightness and/or lack of CO in the lower metaUicity 
environments found at large galactocentric distances. Two SF-idle GMCs with n~10^ cm""^ 
and potentially very low kinetic temperatures (~5K) have been found in M31 (Allen & 
Lequeux 1993; Loinard et al. 1995; Allen ct al. 1995). For Kvir~l-2 (i.e. self-gravitating or 
nearly so) such clouds will have T^q\_q~1.4K and a corresponding Xco~19Xi. In the global 
CO SLEDs of LIRGs such a cold SF-idle component, even if massive, it will be completely 
inconspicious, outshined by the much more CO-luminous SF molecular gas where typical 
t|^o|i_o>15 K. The global Xco of a LIRG, computed using models of its global CO SLED, 
will then invariably be biased by the average ISM properties of its central starburst, often 
having a genuinely low X^o^~(l/3-l/5)Xco,Gai factor (e.g. Papadopoulos & Seaquist 1999; 
Papadopoulos & Allen 2000). Thus its adoption for the entire LIRG can much underestimate 
the contribution from a massive cold, SF-quiescent component where a Galactic Xco applies, 
as indicated by submm, HI and CO J=l-0 imaging for such galaxies. 

We can compute an indicative X^o"^''^ in such a LIRG for a given pco"*^^ by setting 
xS=Xco,Gai for the cold, extended, SF-quiescent component, and X^o^=(l/5)XGai,co for the 
warm, nuclear, star-forming one. Using the results by Papadopoulos & Seaquist (1998) 
where "stacked" CO(2-1)/(1-0) ratios for LIRGs were used to determine a generic warm 
versus cold molecular gas distribution in their disks (their Figure 6) we set Lh~0.9kpc as 
the diameter of the warm SF component and Ls~20 kpc as that of the entire CO-bright gas 
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distribution in the "average" LIRG. With a warm/cold CO brightness temperature ratio of 
tb=T(J!i-o/TS,i_o~3itis: p^lr^^=(l/tb) x [(L,/Lh)'-l]~164. Thus x(rP'^)~xS=X,o,Gai, i.e. 
dominated by the cold disk component. Given that pco'^^ can very greatly among LIRGs, it 
is instructive to also obtain xio "^^^ in a LIRG where the relative distributions of the warm SF 
versus the cold SF-idle gas are well known. For the Sy2 NGC 1068: Lh~2.7kpc and Ls~6kpc 
(Papadopoulos & Seaquist 1999), while tb~5. Thus pco '^•'~0.79 and Xco ^'^''~0.6Xco,Gab 
which is ~3x higher than Xco"* (that would be deduced from its global CO SLED). 



3.4. Moleculeir gas inventories in LIRGs: critical observations 

Highly turbulent molecular gas in the merger-driven starbursts of ULIRGs, and ex- 
tended, cold, SF-quiescent gas in the disks of isolated or slightly perturbed LIRGs represent 
the two extremes where current CO line studies may systematically underestimate the total 
molecular gas mass. In both cases global CO SLEDs often show a low-excitation component, 
seen as small CO J=l-0, 2-1 hne flux "excess" (often with a subthermal CO (2-1)/ (1-0) 
ratio) on top of the dominant CO line emission of a warm and dense SF gas phase with high 
CO hne excitation to much higher J levels. In ULIRGs such a diffuse low-excitation phase 
is typically unbound, warm, and contains little mass. In the disks of LIRGs it consists of 
SF-quiescent, gravitationally bound (or nearly so). Galaxy-type GMCs extending beyond a 
central starburst. In the latter case spatially resolving the CO SLEDs and/or dust emission 
of the nuclear SF region versus the extended cold disk can yield proper molecular gas mass 
estimates not dominated by the (usually) low Xco of the CO-bright nuclear SF region. In 
practice CO and ■'^^CO J=l-0, and J=2-l imaging of those SF and SF-idle areas of LIRGs is 
adequate to determine the nature of the residing ISM (SF: r2i~l, Rio,2i~10-15, SF-quiescent: 
r2i~0.4-0.6, Rio,2i~3-6, the corresponding Xco values, and the areas over which they apply. 

In the upcoming era of ALMA such spatial separations of SF and non-SF ISM in LIRG 
disks along with estimates of the corresponding molecular line ratios and Xco values will be 
straightfoward. For ULIRGs though the low-density and the high density gas components 
can be concomitant or very closely associated, especially if the former is the outcome of the 
disruption of CMC outer layers. Thus high resolution CO observations, even with ALMA, 
may be unable to separate their distributions and their corresponding CO line ratios, leaving 
their Xco factors still uncertain. In such galaxies, with the bulk of the molecular gas mass 
at n>10^cm~^, observations of rotational transitions of heavy rotor molecules such as HCN 
are of paramount importance, even without spatial information. For example a high global 
HCN/CO J=l-0 brightness temperature ratio of ~l/4-l/6 would immediately indicate an 
unsually high Mdense/Mtot(H2), contrasting the ~10x lower such ratio in disk GMCs. Then 
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multi-J observations of heavy rotor molecules can be used to determine the corresponding 
Xmob Xco^ factors and eventually Mtot(H2). 

High resolution interferometric imaging of at least two rotational transitions of a heavy 
rotor molecule (e.g. HCN J=l-0, 3-2) along with at least one of its isotopologues (e.g. 
H^^CN J=l-0) remains an invaluable resource for better determining the distribution and 
mass of the dense gas in the compact disks of ULIRGs. This can be achieved using radiative 
transfer models of the emergent HCN line emission as a function of position within these 
disks, as done earlier for the CO J=l-0, J=2-l interferometric study by DS98. The focus will 
now be on the dense gas where much of the molecular gas mass resides, while its all-important 
average dynamical state (i.e. the Kvir) will be determined as part of the modeling. Imaging of 
the rare isotopologues can much reduce the radiative transfer modeling degeneracies affecting 
Kvii.(r), and yield improved constraints on the dense gas surface densities Sdense(r) of the 
gas disks in ULIRGs. Their total gas mass can then be determined from integrating the 
resulting Edense(i") over the best-fit disk models. 

Here we must note that any molecular hue observations that involve the [C/^'^C] iso- 
tope ratio such as CO/^^CO or HCN/H^^CN, while necessary for reducing LVG modeling 
degeneracies and better constraining Kyir, they involve the additional assumption of the 
[C/^^C] abundance ratio (which we assumed to be 50 in our LVG models). The latter can 
be particularly uncertain in ULIRGs where their large CO/^^CO line ratios have also been 
attributed to a higher [C/^^C] abundance than in spiral disks, a result of early starburst 
ages and/or accretion of relatively unprocessed (by star-formation) molecular gas (Hcnkel & 
Mauersberger 1993). If enhanced [C/^'^C] abundances are indeed the norm in ULIRGs this 
will have the general effect of increasing the various Xco, Xhcn factors since the deduced 
Kvir values per gas phase can now be lower. This will be so simply because larger [C/-'^^C] 
abundances rather than low CO line optical depths (and thus high Kvi,., see Equation 11 in 
Paper 1) can be also responsible for the observed high CO/^^CO hue ratios in ULIRGs. 



3.4-1- Molecular gas mass estimates of ULIRGs: the promise of Herschel 

In ULIRGs a massive, warm, and dense gas phase can have a luminous CO SLED that 
remains prominent up to very high-J rotational transitions. Since J=l-0, 2-1 transitions 
can have significant contributions from a low-excitation diffuse gas component containing 
small fraction of the total molecular gas mass, this leaves only J=3-2 and higher-J CO lines 
as useful probes of the bulk of the molecular gas mass and its properties. The SPIRE/FTS 
aboard the HSO can provide access to CO lines from J=4-3 up to J=13-12 for local luminous 
ULIRGs (e.g. van der Werf et al. 2010; Rangwala et al. 2011) and thus yield a critical dataset 
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for obtaining better total molecular gas mass estimates in such systems. 

Once a fully-sampled global CO SLED from J=l-0 up to J=13-12 is available, a multi- 
component analysis can decompose it into a series of gas components, constraining their 
properties, and using them to obtain the corresponding Xco factors. Such multi-phase mod- 
els could eventually be produced by theoretical advancements in GMC-sized numerical sim- 
ulations (see 3.1.1). Thus HSO observations of high-J CO lines will allow better inventories 
of molecular gas mass in ULIRGs since much of that mass resides in a dense, and presum- 
ably SF and warm phase. We note however that serious degeneracies will remain, especially 
regarding the all-important dynamical state of each gas component upon which the corre- 
sponding Xco factor strongly depends. Multi-J observations of the much fainter (especially in 
ULIRGs) ^^CO isotopologue lines, necessarily using the much larger mm/submm telescopes 
available from the ground (and thus limited by the atmosphere up to J=6-5, 7-6) are im- 
portant for reducing those degeneracies, and better constraining Kyir, or equivalently (for a 
given density and temperature), the escape probabilities per gas component that enter 

the expressions of the Xco factors (Equations 3 and 6). 

4. Probing the extremes: individual (U)LIRGs, their molecular gas, and Xco 

Several (U)LIRGs in our sample merit an individual study either because a larger than 
average number of available CO lines permits it (see Table 7 in Paper I), and/or because very 
high line excitation make their CO SLEDs irreducible to superpositions of star-forming and 
(non-SF) molecular gas (sec also Paper I for a discussion on the implication about ISM power 
sources). Finally there are LlRGs whose global CO(J-|-1-J)/(1-0) ratios suggest the cold SF- 
quiescent clouds found at large galactocentric distances in the Galaxy and the disk of M31. 
Such globally "cold" yet star-forming galaxies are very few, as expected for an IR-selected 
(and thus SFR-selected) LIRG sample, and thus also deserve a closer look as they represent 
the low-excitation range in LIRGs. The detailed model(s) per galaxy and the associated 
discussion are in the Appendix, while here we summarize the most important findings. 

The extreme range of CO SLEDs found for the molecular gas reservoirs of LIRGs dis- 
cussed in Paper I is now marked by individual objects. On the high end are galaxies like 
IRAS 00057+4021, Arp299, IRAS 12112+305 and others whose extremely high CO hue ex- 
citation implies large amounts of very dense (~(l-3)x (10^-10^) cm~^ and (often) very warm 
(Tkin>100K) gas. The large fd=M(n>10'' cm-3)/Mtot(H2) deduced for such (U)LIRGs in- 
dependently recovers, using mid/high-J CO lines, a well-known result for merger-driven 
starbursts obtained using the HCN/CO J=l-0 ratio as an fd proxy (Gao & Solomon 2004). 
These earlier studies (see also Gracia-Carpio et al. 2008; Krips et al. 2008) found ~10 
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times higher dense gas mass fractions in ULIRGs than in isolated spirals. Isofar as the dense 
molecular gas phase closely tracks the dense and warm gas associated with SF sites, the 
mid-J/high-J CO lines of ULIRGs are also expected to show clear indications of a high fd, as 
it is indeed found, and quite unlike what is expected from typical GMCs (fd,GMC~0.02-0.03). 

Furthermore, the very warm and strongly unbound (Kvir>20) states uncovered for the 

massive dense gas reservoirs of some LIRGs can make their global CO SLEDs surpass even 
those expected for SF "hot" -spots in the Orion A and B clouds. While in some distinct 
cases this can be due to strong AGN feedback (e.g. Mrk231), for all other galaxies with such 
extreme ISM conditions (e.g. IRAS 00057+4021, IRAS 08572+3915, IRAS 23365+3604) the 
cause is unclear. It remains to be explored whether the much higher SFR densities of the 
compact SF regions of ULIRGs (e.g. Sakamoto et al. 2008) can create such extraordinary 
ISM states where turbulence-injection by SNRs no longer remains confined in small regions 
but encompasses much more molecular gas mass (see discussion in Paper I). This, along with 
CRs (also SNR-generated) may set up powerful global mechanisms that can volumetrically 
heat large amounts of gas, unhindered by the large dust extinctions and high average gas 
densities that will keep PDRs very localized around SF sites. 

On the low excitation end only a few cases of low global CO line ratios are found 
(e.g. IRAS 05189-2524, IRAS 03359+1523). This is expected for IR-selected (and thus 
SFR-sclcctcd) galaxies (sec Dunne et al. 2000 on the limitations of such samples in finding 
(cold-ISM)-dominatcd systems). In Arp 193 a low-excitation ISM state and the lack of large 
amounts of dense and warm gas is actually suggested by HCN rather than CO lines. Indeed 
despite the presence of a young merger with substantial IR luminosity and a significant 
HCN/CO J=l-0 hue ratio, its global HCN hue emission is consistent with the absence 
of a dense gas phase, possibly a case of strong SF feedback momentarily dispersing its 
dense gas supply (Papadopoulos 2007). Such IR-luminous/(dense-gas)-deficient galaxies will 
be rare (see discussion in B.17) and thus valuable for studying the effects of SF feedback 
onto the dense ISM where the initial conditions of star formation are set. The case of 
IRAS 05189-2524 on the other hand stands out as one where a massive cold molecular 
gas-rich disk is implied but unlike other LIRGs (e.g. IZwl) there are no morphological 
indications whatsover for the presence of such a disk in this very compact ULIRG/AGN 
system. Moreover its "warm" CO(6-5)/(3-2) ratio, warm IR "colors", and compact size in 
cm, near-IR, and optical wavelengths would argue for the presence of only a warm, dense, SF 
gas phase, and thus this object represents a cautionary tale about such conclusions drawn for 
similar high-z ULIRG/AGN systems using only high- J CO lines (e.g. Tacconi et al. 2006). 

We note that massive cold molecular gas disks are imphed for other LIRGs as well (e.g. 
IZwl, VII Zw 031, NGC7469) but without strongly affecting their global CO SLEDs or 
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CO/^^CO line ratios which remains dominated by the warm SF phase. This simply mirrors, 
for molecular lines, a result well-known for dust continuum emission, though in practice 
global CO SLEDs arc somewhat more sensitive to the presence of cold diffuse gas than the 
global dust SEDs are to the concomitant cold dust mass (Papadopoulos & Allen 2000). 

4.1. Effects of average ISM conditions on the Xco 

In Table 1 we tabulate the total molecular gas mass estimates as produced by Equations 

3, 10, the corresponding Xco factors, the minimum molecular gas mass implied for Eddington- 
limited star formation (section 2.3), and the mass of the high-excitation (h)-phase when a 2- 
component model is used to interpret the CO, ^'^CO lines (sec Appendix). A mere inspection 
of the third column demonstrates than in several galaxies (~40%) significantly larger Xco 
factors than the so-called (U)LIRG values of ^^(0.6-1) Xi may apply. For ULlRGs this is 
due to large and even dominant fractions of their molecular gas mass being at much higher 
densities than in disk-dominated LIRGs while for the latter because extended cold molecular 
disks (with a Galactic Xco) can contain much of the molecular gas mass while remaining 
incospicious in the global CO SLEDs used to constrain the global Xco. For ULIRGs this 
low-mass bias includes some very well-known systems such as Mrk231 and Arp220. For 
lower IR luminosity galaxies underestimates of their total molecular gas by the Xco deduced 
from one-phase radiative transfer models of their CO SLEDs can be important in disk- 
dominated systems (e.g. NGC7469), but also in seemingly compact star-forming galaxies 
that nevertheless have "cold" CO SLEDs and Galactic Xco values (IRAS 05189-2524). 

The computed Xco"^'^'' and M^q^^'^'' in Table 1 also make clear that in the case of ULIRGs 
global molecular line SLEDs can reveal the aforementioned mass bias, provided that they 
include mid/high-J CO and/or heavy rotor molecular lines (e.g. HCN). For less IR-luminous 
galaxies with nuclear starbursts and cold SF-quiescent gas-rich disks global SLEDs cannot 
easily identify the presence of the latter except in a few cases (e.g. IRAS 05189-2524). 
This becomes possible only when additional spatial information (e.g. CO 1-0, submm dust 
continuum, or cm imaging) is available for such a disk (e.g. NGC7469). In the absence of 
such information even good one-phase models of global CO SLEDs for a (starburst)-l-(cold 
disk) system may be unable to reveal the disk component leaving the deduced average Xco 
dominated by the starburst phase (e.g. IRAS 02483+4302). In the era of ALMA routine 
CO muIti-J and ^^CO line imaging will provide much better molecular gas mass estimates 
for disk-dominated LIRGs with strong ISM excitation gradients (see 3.4). For ULIRGs on 
the other hand multi-J observations of heavy rotor molecules (e.g. HCN, CS) and their 
rare isotopologues are necessary for confirming the high total Xco and molecular gas masses 
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implied by our results in Table 1, by constraining the all-important dynamic state of the 
massive dense gas phase as discused in 3.2.1. 

5. Conclusions 

In this work our large CO, ^^CO line survey of Luminous Infrared Galaxies (LIRGs) 
detailed in Paper I (Papadopoulos et al. 2011) is used to examine the impact of the wide 
range of average ISM conditions found for these galaxies on their total molecular gas mass 
estimates via the so-callcd Xco=M(H2) /L'qq factor. Our sample includes some of the most 
prominent local ULIRGs (e.g. Arp220, Mrk231, IRAS 17208-0014), often used as templates 
for merger-driven starbursts at high redshifts, as well as less IR luminous disk-dominated 
galaxies. Wc find that one-phase radiative transfer models of the global CO, ^^CO line ratios 
yield (Xco)~(0.6±0.2) (Kkms^^ pc^)^^, similar to that obtained by past studies. The 
average gas temperature and density strongly influence Xco, but the gas average dynamical 
state is the most important influencing factor with unbound gas corresponding to low 
values while self-gravitating gas to higher ones. 

Nevertheless higher Xco~(2-6) (Kkms~^ pc^)^^ values arc deduced for (U)LIRGs 
whenever adequate molecular line data exist to determine the mass contribution of gas at 
densities n>10^cm~^ (high- J CO lines from our survey and/or HCN lines from the litera- 
ture) . Theoretical expectations for the highly turbulent molecular gas in the merger-induced 
starbursts of ULIRGs indicate that, with most of the gas at such high densities, the afore- 
mentioned large Xco values maybe the rule rather than the exception in such systems. Past 
observational studies were unable to determine this, yielding instead much lower Xco values 
(and often considered as ULIRG-appropriate standard ones), because the molecular lines 
used (mostly CO J=l-0, 2-1) could not constrain the properties of the dominant (i.e. the 
dense) gas phase in ULIRGs. Our results indicate that only high-J CO lines and multi-J 
observations of heavy rotor molecules (e.g. HCN, CS) can overcome this mass bias, plac- 
ing the Herschel Space Observatory and ALMA front and center in the quest for improved 
molecular gas mass estimates in ULIRGs in the local and the distant Universe. Of partic- 
ular importance are good constraints on the dynamical state of the dense molecular gas in 
such systems (self-gravitating or not?) since strongly unbound states for this phase seem to 
be the only possible way that Xco in ULIRGs could be much lower than a Galactic value. 
Observations of high-density tracing rare isotopologucs (e.g. high-J ^^CO, H^^CN) will be 
of crucial importance in yielding such constraints. On the theoretical front, CMC-sized nu- 
merical simulations of the turbulent molecular gas with ULIRG-type of boundary conditions 
are needed in order to obtain the dense gas mass fraction and its temporal evolution. This 



-26- 



in turn can inform galaxy-sized simulations of merger systems which recently have included 
molecular gas, but cannot track the dynamical state and mass of the dense gas. 

We also find LIRGs where underestimates of their total molecular gas are the result of a 
massive, cold, gas-rich, disk existing beyond their central starburst. In such cases the global 
CO SLED, and the Xco factor determined from it, remains dominated by the central starburst 
and its often low Xco- Spatially resolving the CO, ^^CO line and dust continuum emission of 
such disks in LIRGs with strong ISM excitation gradients is indispensible in accounting for 
their total molecular gas mass. In such cases, provided that adequate resolution is employed 
to separate the cold gas disk from the nuclear starburst, even low-J CO and ^'^CO lines 
(J=l-0, 2-1, 3-2) are adequate to yield much improved molecular gas mass estimates. 

Finally our study unfolds the wide range of the global CO SLEDs of LIRGs presented in 
Paper I over a subsample of individual systems. We find galaxies whose extreme CO SLEDs 
indicate ISM conditions that surpass those expected for star-forming regions and suggest 
ISM energy sources other than photons from PDRs (as already discussed in Paper I) . These 
can be AGN, extreme turbulence, and/or very large cosmic ray energy densities. Moreover 
we find LIRGs where extreme SF feedback may have momentarily fully disrupted their dense 
molecular gas reservoirs, and LIRGs where large amounts of cold SF-quiescent molecular gas 
are present despite the absence of an extended disk, their compact near-IR/cm emission 
size, and "warm" CO (6-5)/(3-2) ratios. The latter type of objects can be particularly 
worrying if encountered at high redshifts where such characteristics can readily lead towards 
large underestimates of their total molecular gas mass. The wide range of the average 
ISM conditions and the intriguing possibilities that may lie behind it, make our subsample 
of individually-studied LIRGs an excellent target for future ALMA and NOEMA (for the 
northern objects) molecular line imaging observations leading towards a complete picture 
about ISM energetics, AGN and SF feedback on the molecular gas in galaxies. 
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night inspirations, and to his wife Map'-fapira for her enduring support. The project was 
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A. The Xco factor in an LVG setting 

The CO J=l-0 line luminosity 'L'qq(1-Q) (in Kkms~^pc^, see Equation 5, Paper I) can 
be re-expressed as 
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/ / Tb,i_o(F,V)dadV= / / 
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where any given velocity V is assumed to correspond uniquely to a source surface As(V), 
emitting at Tb(f, V). In turn this can be parametrized by a "depth- in-the-source" parameter 
R(V) so that these iso- velocity surfaces completely scan the entire source volume for a range 
ARs (which corresponds to the FWZI of the source velocity field) , and without any radiative 
coupling (the LVG assumption). The latter implies that the line luminosities emanating 
from these surfaces simply add up. Thus we can write 
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where (..) denotes averaging over the entire source volume AVg, and for an assumed constant 
velocity gradient. Thus for 
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and M(H2)=/xm(H2)(n(H2))AVs {/J,— 1.36 accounts for the He mass) Equations Al, A2, and 
A3 along with subsituting astrophysical units yield for the Xco factor 
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The last expression (used in the main text with the averaging symbols omitted for simplicity) 
for a=1.5 becomes identical to that derived by Solomon et al. 1997 (their Equation 21 for 
f=l). For virialized gas motions (Kvii.=l), and typical conditions in Galactic GMCs with 
(n(H2)) = (100-500) cm-3 and Tb,i_o=10 K yields Xco~(3-6) M0(Kkms-i pc^)"S in good ac- 
cord with the average Galactic value and its uncertainties. Finally, expression A4 is valid for 
any optically thick molecular line emission used as a mass tracer of a particular gas phase 
(e.g. HCN J=l-0 tracing dense gas with n(H2)>10^cm~"^), while multi-J observations of the 
particular molecule and its isotopologues can be used to constrain (n(H2)), Tb,i_o, and Kyir. 
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A.l. The optically thin approximation for Xco 

From standard formalism the integrated line luminosity for an optically thin CO J+l^J 
line, omitting the CMB for simplicity, is given by 

Lj+i,j = Nj+i Aj+i jhi/j+i,j, (A5) 

where Nj+i is the total number of CO molecules at the J+1 state and Aj+u is the Einstein 
coefficient of the J+1— >J transition. This line luminosity can be re-expressed in terms of 
Lj^^j (in Li^Kkms"-*^ pc^ units, see Al) using 

87rkjj7/''J J , 

Lj+i,j = Lj+i,j, (A6) 

which combined with A 5 yields 

87rkB^j\i,j , 

hc3Aa+,/j+^J- ^^^^ 
The total molecular gas mass is then given by 



M(H2) = (No + Ni + N2 + N3 + N4....)Rco/iniH2, (A8) 

where Rco=[H2/CO] = 10'' is the CO abundance and /x=1.36 accounts for He mass. The 
molecule population at the J=0 level obviously cannot be estimated from a line transition, 
and we thus compute it from No=(go/gi)exp[Ei/(kBTex,io)]Ni where gj=2J + 1 denotes the 
J-state degeneracy factor, Ei/kB~5.5K, and Tex,io is the excitation temperature of the CO 
J=l-0 transition. Substituting this in A8 yields 

M(H2) = (^1 + le^-^/^-- + ^ + ^ + ....^ Rco/iniH.Ni, (A9) 
which after substituting the expressions from A7 yields 
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where rj+u are the Lco(J+1-J)/L'co(1-0) ratios observed (Equations 5, 6). The Einstein 
coefficients scale as a function of J as 



Ai-i-i T — 3- — Ai n 
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while i/j+ij=(J + l)uio. Thus equation AlO finally becomes 
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Substituting the various physical constants and introducing astrophysical units yields, 
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In the literature the LTE approximation is often used by setting X]k=o ^k=ZLTE/giNi and 
ZLTE~2(kBTk/Ei)=2[Tk/(5.5K)] in Equation A8, which then yields 
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For optically thick CO emission where the optical depths arise in small gas cells with 
respect to molecular cloud sizes CO line emission remains effectively optically thin (i.e. 
traces the entire emitting gas mass) but with Aik— >AkAik in all the previous equations, with 
Ak=[l — exp(— /3ik)] /Tik being the photon escape probability (for a spherical cloud). In such 
a case Equation AlO can be trivialy modified for finite CO line optical depths to yield 
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B. Models of individual LIRGs 

Here we describe the radiative transfer models in selected galaxies from our sample. 
The details of the Large Velocity Gradient (LVG) radiative transfer code used, and the 
parameter space explored by its [Tkin, n, K^ir] variables can be found in Paper I. In the several 
cases where a one-phase gas LVG model manifestly fails to reproduce the available line 
ratios we use a two-phase approximation (see 2.4) to better represent the underlying average 
molecular gas properties, and compute the corresponding Xco"^*^^ (Equation 10). The (h)- 
phase properties and its Xco^ are constrained using CO J=4-3, 3-2 and ^^CO lines, with 
the 6-5 transition placing constraints as a lower limit of the true emergent line luminosity 
(see Paper I), with higher CO J=6-5 luminosities typically yielding even denser gas LVG 
solutions with larger Xco"*- The "residual" CO line ratios (determined from Equations 8, 9) 
are used as inputs into our LVG radiative transfer code in order to constrain the physical 
properties and corresponding X^^q factor of the low-excitation (l)-phase. In the cases where 
a two-phase ISM model is necessary but inadequate CO line observations exist to constrain 
the (h)-phase, we assume one based on the model defined by the Orion star-formation "hot- 
spots", their SLED, and molecular gas mass normalization (see 2.4). Finally in the very 
few cases where heavy rotor molecular lines (mostly HCN transitions) exist, they are often 
used to determine the (h)-phase. For all the LIRGs studied, their measured CO line ratios, 
computed SF-powered IR luminosities Ljj^ , and the dust temperatures used to constrain the 
LVG model (i.e. Tkin/Tdust^l) can be found in Paper I. 



B.l. IRAS 00057+4021 

There are only a few studies of this LIRG, and the only available CO J=l-0 interferomet- 
ric map shows a very compact nearly face-on source with most of the gas in a < 1" (~864pc) 
core (DS98). Its luminous CO J=4-3 line emission cannot be accounted by one-phase LVG 
models, with the best such model: Tkin=90 K, n(Il2)~3 x 10^ cm~^, Kvir=l, reproducing well 
the observed T21, and the R21 lower limit, but yielding r43'^''=0.48, which is ~l/3 of the 
observed value. Moreover the fact that r43>r43'' (Orion) (sec 2.4) and r43>r32, makes obvi- 
ous that a superposition of low-excitation (SF-quiescent) and high-excitation Orion-type SF 
"hot" spots cannot produce such a highly-excited global CO SLED, and that even extreme 
GMCs consisting solely of such "hot" spots are inadequate. Thus the average state of the 
molecular gas in IRAS 00057+4021 "out-excites" even the SF regions in Orion A and B. 
The existence of such ULIRGs, and the implications for an ISM energy input other than 
far-UV/optical/IR photons from SF sites (i.e. turbulent and/or CR heating) has already 
been noted in Paper I (see also Rangwala et al. 2011). 
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Using 143-17=0.89 along with the other CO hne ratios still does not yield an average ISM 
state that can reproduce them. The best LVG solution ranges then arc: Tkin=(35-65) K, 
n(H2)=103cm-3, Kvi,~l-4, and Tki„=(70-150) K, n(H2)=3 x lO^cm-^, K^i^-l. These re- 
produce the observed r^2, r2i and the R21 lower limit, but yield r43(lvg)~0.45-0.50, i.e. still 
significantly lower than i^s-cr. A 2-phasc model using an Orion-defined (h)-phase (and SLED) 
gives r2i = 0.91, rS=0.53, r2=0.83 (Equations 8, 9), and remains unable to fit these residual 
hne ratios, failing to account for the J=4-3 line (r43(lvg)~0.3-0.4). Using only the J=3-2, 
J=4-3 transitions as constraints (r43/32 =2.23 ± 0.97) of the (h) phase and restricting LVG 
solutions with Tkin>Tdust(=37K) and Kvir>l, yields Tkin>95K, n(H2)=3 x 10^ cm^^ and 
Kvir~20. The fit yields no upper hmit on Tkm, as expected since the maximum possible r43/32 
(attained for optically thin thermalized CO lines) is (4/3)2c~'^:^2/(kBTkin) = i,77c-i6.fi/Tkin_^l,77 
for 16.6/Tkm— >0. Using the lower r43/32=1.26 value (=(obscrvcd)— a) as input in the LVG 
code yields identical results. Thus a very warm, dense, and strongly kinematically stirred gas 
phase is present in IRAS 00057+4021. An estimate of its mass can be obtained from the mea- 
sured Lco(4_3) = (5.84 X 109)Li, and l[5'q(''^''2o) =V(rS~''''VLco(4-3)' where (ri5"''''^)=3.4 is 
the average CO (4-3) /(l-O) ratio computed from the range of LVG solutions compatible with 
the observed one. Thus for L^^Q^^}_f^^ = {1.72 x 10^) Li, and X^o^~0.92Xi (computed from the 
LVG solutions compatible with r43/32) it is Mh_ex(H2)~1.6xlO^ Mq, comparable to the total 
(HI-I-H2) gas mass of the Milky Way and the minimum MsF(H2)=L[^^/eg,*~1.2xlO^ M© ex- 
pected for "fueling" star formation in this LIRG. In Table 1 wc comprehensively list all these 
gas mass estimates for IRAS 00057+4021, along with Mtotai(H2) computed using 1-phase and 
2-phase LVG models (Eqs 3 and 10). We must note however that for IRAS 00057+4021 and 
most (but not all) LIRGs where 2-phase models are used, the large degeneracy that exists for 
LVG solutions of the "residual" CO line ratios of the (l)-phase (cool and near-virial versus 
warm and strongly gravitationally unbound states) translates to a large range of xS, and 
thus X^o"^''^ values (Table 1). This in effect is the low- J CO SLED degeracy described in 
3.4, which the global CO SLEDs available in our study cannot "break" . 

The high-excitation (h)-phase in IRAS 00057+4021 contains ~(20-50)% of Mtot(H2) 
(Table 1), which is at least ~10x larger than what would be expected if its ISM was reducible 
to typical Galactic GMCs. Moreover its r^i^(lvg)~2.5, R^i^(lvg)~23 ratios (for Tkin=100K, 
n(H2)=3xl0^cm~^, and Kvir=22) are higher than those in the Orion SF "hot-spots" near 
HII regions and 0,B star associations (Sakamoto et al. 1994). This further corroborates the 
irreducibility of the observed CO SLED of IRAS 00057+4021 to ordinary photon-powered 
ISM states. One the other hand strong mechanical feedback from SNR-drivcn shocks could 
in principle yield such highly excited CO SLEDs (Arikawa ct al. 1999; Bolatto ct al. 2003), 
but just like radiative feedback from O, B stars, such SNR-driven "hot" spots involve only 
small fractions of a typical CMC mass. Indicatively, in the Galactic SNRs W44 and IC443 
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a warm, dense, and very kinematically stirred component with r2i~1.3-4 (Seta et al. 1998) 
involves only ~1% of the impacted GMCs, with negligible effect on their global CO line 
ratios ((r2i)GMC~0.6). We must also note however that IRAS 00057+4021 is a Seyfert2 and 
its AGN, if X-ray luminous, can induce such extreme conditions to a large molecular gas 
mass, if distributed close to it (Schleicher et al. 2010). 

B.2. IRAS 00509+1225 (IZwl) 

This LIRG hosts an optically and X-ray luminous QSO in the center of a molecular 
gas-rich disk, the site of a vigorous circumnuclear starburst (Barvainis et al. 1989; Eckart et 
al. 1994 and references therein). The latter has been imaged interferometrically (Schinnerer 
et al. 1998; Staguhn et al. 2004) revealing a disk of ~8-12 (~12-14kpc) with cold ISM but 
increasing ISM excitation towards its central ~1.5 (~2kpc) SF circumnuclear ring. 

In terms of global CO line excitation while its r2i=0.84 ± 0.24 ratio can be (barely) 
compatible with that of a SF-quiescent disk like the Milky Way (~0.6), its luminous J=3-2 
line with r32=1.16 ± 0.40 is clearly dominated by much higher ISM excitation, typical of 
a SF phase (even r32-o"=0.76 is 2.5x higher than r32=0.3 typical for SF-quiescent ISM). A 
one-phase LVG model yields an acceptable fit within the measurement uncertainties, with 
r2i(lvg)~l, r32(lvg)~0.90-0.95, and Rio~8, R2i~5-7 while Kvir~l-2 (a result of the modest 
Rio and R21). The best one-phase LVG solutions yield Xco~1.5Xi (Equation 3), and a 
total molecular gas mass M(H2)~8.5xlO^ M©. The minimum dense gas mass fueling an 
Eddington-limited SF is Msf(H2)=4x 10^ Mq, which only ~5% of the Mtot(H2) and typical 
for ordinary SF GMCs found in galactic disks. Thus I Zw 1 stands as an example of a LIRG 
with star formation occuring in an isolated gas-rich disk, fueled by ordinary molecular clouds. 

The reasonable 1-phase LVG fit of the observed global CO line ratios in I Zw 1 make a 
2-phase decomposition unnecessary. Nevertheless the available observational information on 
the presence of a cold molecular gas-rich disk and a circumnuclear ring starburst allow such a 
decomposition as a useful test on the limitations of global CO SLEDs and their fit by a single 
average ISM state, to yield reliable total molecular gas mass estimates. Using the CO J=3-2, 
6-5 and ^^CO J=l-0 lines to constrain the (h)-phase we find the best LVG solution range 
at Tkin=(65-75) K, n^lO^cm"^ and Kvir=13 that yield X^o^=0.65Xi. Adopting a Galactic 
Xco=5 Xi for the computed residual CO J=l-0 emission, which we attribute to the extended 
gas disk, yields Xco ^^^^=1.65X1 (Equation 10). similar to that computed from the one-phase 
LVG model. We caution however that the CO J=6-5 line luminosity is highly uncertain and 
along with it the current 2-phase decomposition. 
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B.3. 



NGC 828 



This is a disturbed spiral galaxy with a prominent dust lane and Ha emission in its 
center and two sources symmetricaly around the center source along the major axis (Hattori 
et al. 2004). VLA imaging of its cm continuum emission show it elongated along the disk, 
as would be expected for star formation, with no signs of AGN activity such as a radio core 
or jets (Parma et al. 1986). Its single dish CO J=l-0 detection showed abundant molecular 
gas (Sanders et al. 1986) while subsequent interferometric CO 1-0 imaging with OVRO 
recovered all the single dish flux and showed it to be extended ~10"x20" (~7.5kpcx3.6kpc) 
with the longest dimension lying along the optical disk (Wang et al. 1991). The same study 
obtained a dynamical mass of Mdyn~4xl0^° Mq within a CO disk radius of 3.9 kpc (for the 
adopted cosmology). However at ~4" (1.4 kpc) south-east of its nucleus the CO-derived 
velocity field deviates from that of a normal rotating spiral, possibly indicating the presence 
of another molecular gas concentration brought in as part of a merger (Wang et al. 1991). 

The available CO lines are consistent with a star-forming disk with a well-excited CO 
J=3-2 line (r32=0.70±0.15, Paper!) but whose modest Rio=10±2 and R2i = 12±3 ratios indi- 
cate that any merger activity has not disturbed the average dynamical state of the molecular 
clouds. This is reflected in the one-phase LVG solutions found for Tkin>35K which invari- 
ably have Kvir~2-4. These solutions remain highly degenerate however with Tkin=(35-50) K, 
n=10^cm~^, and Kvir~4 almost as good as Tkin=(55-150) K, n=3xl0^cm~^, and Kvir~2. 
The corresponding Xco~(0.8-l.l)Xi, yields Mtot~(4.6-6.3) xlO*^ Mq, while Msf~3x10^ Mq is 
necessary to fuel an Eddington-limited star formation. 

Thus the dense and warm gas fueling SF sites in this low IR luminosity galaxy will 
amount to only ~5%-6.5% of the total molecular gas mass, consistent also with its small 



SLED can still be easily dominated by even small amounts of star forming gas whose of- 
ten modest Xco may be much smaller than that of an incospicious SF-quiescent cold gas 
reservoir. For NCC828 there is apriori knowledge of an extended molecular gas disk with 
~8kpc diameter. However a 2-phase fit with an assumed (h)-phase SLED and mass normal- 
ization (see 2.4) leaves the (l)-phase properties largely undetermined with Xco~(0.7-2.5) Xi 
(~X^o~^''^ since the (l)-phase CO J=l-0 luminosity dominates), and Mtot~(4-14) xlO*^ Mq. 



HCN/CO J=l-0 ratio of r 




0.022 (Gao & Solomon 2004). Nevertheless global CO 



B.4. IRAS 02483-1-4302 



This two- nuclei system is a gas-poor/gas-rich merger of an elliptical (nucleus A) going 
through the disk of a former spiral containing nucleus B (KoUatschny et al. 1992) and all of 
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the molecular gas detected in a CO J=l-0 interferometer map (DS98). The optically bright 
nucleus A hosts an AGN with a Sy2 spectrum, and is located 3.8" (~3.76kpc) to the west 
of nucleus B. The latter is where star formation occurs in this system, "activated" by the 
merger event which also likely triggered the QSO activity in the gas-poor merger progenitor. 

Its CO line ratios up to J=3-2 can be well-fitted by a wide range of conditions with 
even a cold Tkiii=15K but highly unbound (Kvir=22) phase yielding a good fit. Most solu- 
tions are found at warmer temperatures though with Tkin= (35-55) K, and Tkin=(90-120) K, 
low densities of n~(3x 10^-10'^) cm^^, and strongly unbound states Kvii.=7-40. Such a wide 
degeneracy of LVG solutions representing the average conditions of the molecular gas in 
IRAS 02483+4302 is typical for LIRGs where rj+i,j<l, and no CO/^^CO ratios are available 
to constrain hne optical depths (and thus K^ir). For this LIRG a lower hmit of Rio>13 is not 
enough to eliminate the aforementioned degeneracy, though it does preclude SF-quiescent 
and virialized states. For for all the LVG solutions compatible with its global line ratios 
we obtain Xco=(0.25-0.70) Xi (Equation 3). Restricting the admissible range of LVG so- 
lutions with Tkin/Tdust^l, and Kvir<20 yields Xco~0.45Xi and a total molecular gas mass 
M(H2)~1.6x 10^ Mq, similar to that obtained using the optically thin, LTE, approximation 
^thin^''~2x 10^ Mq (Equation 5 for Tkin=55K) a result of the low/moderate optical depths 
of the LVG solutions used to deduce Xco in this LIRG. 

The minimum gas mass for an Eddington-limitcd star formation in IRAS 02483+4302 
is Msf(H2)~1.4x10^ Mq, which is ~90% of its total molecular gas mass (for Xco=0.45Xi). 
This is rather puzzling as the SF phase, with its expected high densities (>10*^cm~^) and 
near gravitationally bound states (Kvir~l-5) can easily have GO lines thermalized and with 
substantial optical depths up to at least J=3-2. Thus if the SF phase indeed dominates 
the total molecular gas mass in this ULIRG, one would expect m|^™^<^M(H2) rather than 
M^^™^~M(H2) as the latter implies small/moderate optical depths for the bulk of the molec- 
ular gas mass. Moreover the physical states compatible with the global CO line ratios of 
IRAS 02483+4302 while highly degenerate, are hardly indicative of a dense SF gas phase 
dominating its total molecular gas mass. It is also telling that all one-phase LVG solu- 
tions fail to reproduce the GO J=6-5 line luminosity of this LIRG, indicating the presence 
of another, potentially massive, gas component whose higher GO line excitation becomes 
prominent beyond J=3-2. Using only its CO J=3-2, 6-5, and the upper hmits on the ^^CO 
J=l-0, 2-1 hues as constraints of the (h)-phase in a 2-phase model still yields a significant 
range of solutions though all have n~(10^-10^) cm~^ with Tkin>60K, typical of a dense and 
warm SF-related phase. Most (but not all) have virial states (Kvir~l) with X^o^~(1.5-2.5) Xi. 
From the corresponding (rg5^)ivg~0.90 and (Xco^)ivg~2.2Xi we obtain L^|J\_q=1.5x10^ Li, 
and Mh-ex=3.3xl0^ Mq. Higher GO J=6-5 line luminosities will only make this mass larger 
(both via a higher l|.q^^_q and a higher deduced (Xco'*)ivg)- 
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The "residual" global CO line ratios of the (l)-phase in this LIRG cannot distinguish 
between a cold, SF-quicsccnt gas at near virial dynamical states (and Galactic Xco) versus a 
warm and highly unbound phase expected in the highly turbulent environments of mergers 
(and with a low Xco because of the high K^ii.). Thus for pco ^^^=1.43 computed from our 
2-phasc model, X^o^=2.2Xi and xS = (0.5-2.5) Xi, we obtain xio"''^^ = (l. 2-2.4) Xi. Higher, 
near-Galactic, values remain possible if high resolution CO imaging were to reveal the (1)- 
phase as an extended cold, SF-quiescent gas reservoir, something that global CO SLEDs 
cannot easily do. 



B.5. IRAS 03359-1-1523 

Two interacting galaxies ~10" (~6.9kpc) apart can be discerned in optical images, with 
only the eastern source being bright in radio wavelengths (Condon et al. 1990; Goldader 
et al. 1997) and containing most (>,75%) of the CO J=3-2 emission (Leech et al. 2010, 
Paper I) . This system has the lowest CO(3-2) / (1-0) ratio (=0.18±0.05) indicating the lowest- 
excitation global CO line excitation in our sample, though we cannot exclude the possibility 
that significant CO J=3-2 flux was "missed" by the 2-point observations of this two-nuclei 
system. 

The best LVG solution indicates SF-quiescent gas (Tkin=15 K n=3xl0^ cm~^ and Kvir~2), 
yielding r32~0.22 and R2i~18 which is somewhat higher than the observed value (=12 ± 3). 
The corresponding Xco~2.5Xi gives Mtot~2.2xl0^° M©, while the LTE approximation for 
Tkin=15K gives M^|j'™''=1.8xl0^ M©, which is ~10x smaller than Mtot as expected given 
the large CO line optical depths (rio>14) of the corresponding LVG solution. The SF-related 
molecular gas mass MgF(H2)=8.2xlO®M0 amounts to only ~4% of M(H2), consistent with 
the low average CO line excitation of this system and ensembles of ordinary GMCs. 



B.6. VII Zw 031 

This galaxy has one of the highest CO J=l-0 luminosities in our sample and while 
early ground-based images were suggestive of an elliptical system (Sanders & Mirabel 1996), 
near-IR NICMOS images clearly show a spiral disk with very bright asymmetric arms tightly 
around its nucleus over scales of R~(l-1.4) kpc (half-light radii), with numerous star clusters 
on its disk (Scoville at al. 2000). This is one of the few LIRGs where there is no evidence 
of what has triggered its starburst activity, namely neither a nearby companion system 
(indicating an early merger stage), nor tidal tails (signs of an advanced or post- merger 
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system). High resolution interferometric imaging indicates a nearly face-on source with a 
rapidly rotating gaseous ring on scales of few hundred parsecs (DS98). 

The low-J CO SLED with r2i=0.72±0.12 and r32=0.43±0.19 indicate rather low aver- 
age hne excitation, compatible even with SF-quiescent ISM (r2i~0.5-0.6, r32~0.25-0.30). A 
radiative transfer model of only the CO J=l-0, 2-1, 3-2, and ^^CO J=l-0, 2-1 hues yields 
Tkin= (30-65) K, n=3xl0^cm^^ and Kvir~2 (near-virial) . The corresponding Xco~1.25Xi 
factor yields M(H2) = 1.5x 10^° Mq while the minimum molecular gas mass needed to fuel 
an Eddington-limited star formation in this LIRG is: Msf(H2)~2x 10^ Mq. On the other 
hand a luminous CO J=4-3 with r43=1.46±0.45 (Paper 1) is highly excited, and the reason 
why one-phase LVG solutions, while reproducing the observed r32,2i and Rio,2i hne ratios 
within their measurement uncertainties, they yield r43'^^<0.20-0.26, i.e. ~6-7 times lower 
than observed. The CO J=6-5 line luminosity with r65=0.22±0.07 is also much higher than 
predicted by these models (r^^g ^^=0.01-0.02). Interestingly the ^^CO J=l-0, 2-1 lines, with 
modest Rio,2i~10 and r2i(^^CO)~0.72 ratios, also indicate a low/average excitation state. 

Thus only the CO J=4-3 and J=6-5 hne luminosities signify the presense of another 
potentially massive warmer and denser gas phase, demonstrating once more the importance 
of high- J CO hne measurements. Moreover, the r43>r32 inequahty sets the global CO SLED 
of VII Zw 031 apart from those reducible to a mixture of a dense and warm gas phase (with 
a CO SLED typical of Orion hot-spots), and a cooler, diffuse one associated with non-SF gas 
and a low-excitation SLED. As in IRAS 00057+4021, this implies large amounts of hot and 
dense gas, maintained in a strongly unbound dynamical state as to keep the average optical 
depth of even high- J CO lines below unity (see also discussion in Paper I). We note however 
that in VII Zw 031 the CO J=3-2 line luminosity is highly uncertain because of the large 
system temperature, though even a 3x stronger CO J=3-2 line would maintain the afore- 
mentioned inequality and its implications for extraordinary ISM conditions in this galaxy. 

We use the CO (4-3)/(3-2) ratio r43/32 to constrain the (h)-phase in a 2-phase model 
while using the CO J=6-5 luminosity only as a lower limit. We also assume the minimum 
i'43/32=(L'co(4-3)"'^)/(^co(3-2)+^)~-'--^'^ valuc consistcut within the measurement uncertain- 
ties in order to place hmits on the lowest possible average excitation level for the (h)-phase. 
This will still be high as even that minimum value is close to the theoretical maximum 
of (4/3)^, attained only for a hot gas phase with fully optically thin and thermalized CO 
lines. We find good solutions only for Tkin>95K which keep improving up to the max- 
imum Tkiii=150K considered, while n=3xl0^cm~^ and Kvii.=22. No good solutions can 
be found for lower temperatures as these have unphysical Kvir<l values. Over the good 
solution range the corresponding (r43'*)=3.4 and rg5^>2. The corresponding (Xco^)=0.9Xi, 
and from L[iQ^_Q=(l/(r43'*))L|-,Q^4_3j~3. 45x10^ Li (we used the C0(4-3)-(T value) we esti- 
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mate: Mh-ex=3.1xl0'^ Mq. As expected, with a ^=^CO(2-1)/(1-0) ratio of ~0.73 (subther- 

mal), an (h)-phase model cannot successfuly include the ^^CO lines , while their modest 
global Rio,2i~10 (Tabic 8, Paper I) are typical of disk GMCs rather than of those in star- 
bursts (where Rio,2i^l5). This along with the fact that there is very little CO J=2-l and 
3-2 emission left after subtraction of the corresponding (h)-phase line emission suggest a cold 
disk-like (l)-phase. Assuming an Xco=4.5Xi (i.e. Galactic), and for pccr'^'*=2.34, we obtain 
X^o~^^^=3.4Xi and a total gas mass of M^^7''^'=3.9xl0-^° Mq. Thus despite its prominence 
in CO J=4-3 and 6-5 lines, the high-excitation component amounts to only ~8% of the 
total molecular gas mass in VII Zw 031. Interestingly, the CO J=2-l interferometric map 
by DS98 resolves out half of the 30-m single dish flux, indicating extended emission beyond 
the ~2"(2kpc) nuclear region where the bright CO J=l-0 emission is distributed. Indeed 
in the DS98 interferometric image faint CO J=l-0 emission can be seen out to a radius of 
~4" (4kpc). New sensitive CO and -"^^CO J=l-0 and 2-1 interferometric observations with 
similar u-v coverage (and thus enabling reliable line ratio maps) can be used to discern the 
existence of such a massive cold molecular gas disk in this intriguing object. 



B.7. IRAS 05189-2524 

This is Sy2 galaxy that appears to be a late stage merger with a very red compact 
nucleus and a tidal tail, (Farrah et al. 2003; Veilleux et al. 2002). This (U)LIRG is the most 
compact and has the warmest mid-IR colors (f25/im/f60/im=0.25) out of a nearly complete 
ULIRG sample selected for their "warm" mid-IR colors (f25/im/f60/im>0.2). Such systems are 
are thought to represent a critical stage in the so-called ULIRG— > (optically luminous QSO) 
transition scenario (Sanders et al. 1988b, c). The AGN in this ULIRG is X-ray luminous 
with Lx(2-10keV)~1.3xl0''^ergss-^ (Dadina 2007). NICMOS near-IR imaging reveals an 
unresolved nucleus in all three near-IR bands, with a half light radii of ~(100-140) pc while 
CO J=l-0 observations show large amounts of molecular gas (Sanders et al. 1991) but no 
interferometric CO images are available yet for this southern ULIRG. 

Unlike all other mergers/starbursts, IRAS 05189-2524 has a small R2i=6 ± 2 ratio, 
typical of the self-gravitating (or nearly so) GMCs in spiral disks. As a result all LVG so- 
lutions within the expected temperature range (i.e. Tkin/Td>l) have Kvir~l-2, indicating 
virialized gas motions, quite unlike the much higher Kyir found in merger/starbusts. The 
observed r2i=0.67 ± 0.15 together with the modest R21 arc actually perfectly compatible 
with a SF-quiescent phase (where r2i~0.6), although a SF-active one (r2i~0.8-l) is certainly 
not excluded within the measurement uncertainties. The best LVG solutions converge to- 
wards conditions typical of Galactic GMCs but at elevated temperatures with Tk=(30-55) K, 
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n(H2)=3x 10^ cm-^ and Kvir~l and yielding r2i=0.86-0.91, 132=0.50-0.67, 143=0.30-0.40 and 
R2i=9. Using values within the rj+i j ± a range as constraints does not change the basic 
picture of low-density, somewhat warm gas, in a virializcd or nearly so dynamical state. For 
the best solution it is: Xeo~3.5Xi, which yields Mtot(H2)=l. 37x10^° Mq. 

Denser and warmer gas must of course be present, associated with the vigorous star for- 
mation in IRAS 05189-2524 and the minimum mass needed to fuel it of Msf(H2)~3.6 x 10^ M© 
The CO J=6-5 line also indicates gas with significantly higher excitation than that in- 
dicated by the lower-J CO lines since, for the optimal LVG solutions for the low-J CO 
SLED segment: r65(lvg)~0.01-0.05, much lower than the observed r65=0.37 ± 0.14. Us- 
ing a 2-phasc model with an Orion-derived SLED and an Eddington-limit mass normaliza- 
tion for the (h)-phase (see 2.4) yields (l)-phase ratios of r^\^=0.45, r32=0.28 and R2i=3.5 
with the best LVG solution being Tkin=15K, n~300cm~^ and Kvir~2.2, i.e. typical of SF- 
quiescent GMCs found in the disk of the Galaxy. The corresponding Xco=2.5Xi, which for 
a computed Lqq^^_q^=2.3x10^ Li, yields a total mass of cold SF-quiescent molecular gas of 

^GxIO^Mq. For the estimated p^o"''^=1.43 and X^o^=2.2Xi (see 2.4) we find X^o"P''^~2.4Xi 
and a M^q7^^^~9.4x10^ Mq with ~65% of this gas mass in a cold SF-quiescent phase. 

This merger ULIRG/AGN with its large molecular gas reservoir, nearly-Galactic Xco 
(~3-6 times higher than the so-called ULIRG- values of ~(0.6-0.8) Xi), and most of its molec- 
ular gas mass in a cold reservoir, stands as the clearest counter-example to standard views 
regarding this class of galaxies in the local Universe. Moreover, unlike I Zw 1 (another promi- 
nent "warm" LIRG considered midway in a supposed ULIRG— > (optical QSO) transition) 
whose known spiral disk can be the repository of its large SF-quiescent molecular gas reser- 
voir, the very compact size of IRAS 05189-2524 in cm continum, near-IR, and optical wave- 
lengths (Condon et al. 1990; Scoville et al. 2000; Surace et al. 1998) seems to preclude such 
a configuration. This compactness along with a "warm" CO(6-5)/(3-2) ratio of ~0.63 (i.e. 
dominated by the SF phase) in a ULIRG/AGN whose molecular gas reservoir is dominated 
by a cold SF-quiescent component offer a cautionary tale regarding similar ULIRG/AGN 
objects at high redshifts where such charecteristics were often used (and still are) to argue 
for the presence of only a dense and warm SF gas reservoir (Tacconi et al. 2006). Future 
ALMA observations of this southern ULIRG and similar objects would thus be very inter- 
esting for revealing the distribution of their large molecular gas reservoir and especially of 
their cold phase. 
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B.8. IRAS 08030+5243 

This LIRG is one of the very few (=3) appearing as a single undistorted nucleus in a 
large near-IR (2.2yum) and optical imaging survey of 56 ULIRGs selected from the IRAS 2 Jy 
sample (Murphy et al. 1996). Its CO line ratios reveal a low average excitation state but lack 
of ^^CO observations permit a very wide range of conditions to reproduce its low global 
and r2i ratios. In turn these conditions yield a large range of Xco values, from Xco~4.6Xi (for 
Tkin=15K, n~300cm-3 and Kvir~l), to Xco~(1.5-1.7) Xj (for Tkin=(55-90) K), n^lOOcm-^ 
and Kvir~l) and also down to Xco~0.61 Xj (for Tkin=[20K, (40-50) K], n^[10^, 300 cm'^] and 
Kvir~[13, 7] respectively). For a set of very warm (Tkin=(95-130) K) and strongly unbound 
(Kvir=22) diffuse gas (n=3xl0^ cm~^) solutions the Xco factor can become as low as ~0.3Xi. 
It is worth remembering that such gas may indeed exists in most ULIRGs as an "envelope" 
phase of a much denser dense phase, a result of disrupted GMCs during a strong merger. Its 
high brightness in low-J CO lines (a result of its high Tkin and Kvir values) is responsible for 
the low Xco deduced for such systems (which nevertheless can have most of their molecular 
gas mass in a much denser phase with high Xco, see discussion in 3.2). 

Demonstrating the importance of ^'^CO line measurements we note that the CO/^^CO 
J=l-0 ratio changes dramatically among the aforementioned LVG solutions groups, from 
Rio~5 for cold virialized gas, to Rio~30 for the warm highly non-virial ISM states. Never- 
thelesss lack of such measurements makes it impossible to constrain the LVG solution groups 
and thus Xco in this ULIRG. Values as low as Xco=0.3Xi are likely excluded however as they 
yield Mtot~3x 10^ Mq, comparable to MsF(H2)=2.4xlOni0. In such a case the global CO 
SLED would be that of dense warm SF gas (i.e. highly excited) rather than of low-excitation. 

B.9. IRAS 08572+3915 

This is another prominent ULIRG from the "warm" ULIRG sample (Sanders et al. 
1988c) consisting of a close pair of interacting spirals with their nuclei separated by 5" 
(~5.6kpc) and clearly discernible in cm, near-IR and optical images (Condon et al. 1990; 
Scoville et al. 2000; Surace et al. 1998). Its NW nucleus is unresolved in near-IR and is 
the only one detected in CO J=l-0 interferometer maps where it appears unresolved with 
^co<2.1"(~2.3kpc) (Evans et al. 2002). Mid-IR (12-25 /xm) imaging showed all the mid-IR 
emission also emanating from the NW nucleus and a region <330 pc in diameter (Soifer et 
al. 2000, and for the adopted cosmology). 

The discovery of very luminous CO J=6-5 emission towards the NW nucleus (Paper I, 
Figure 2), corresponds to second highest tq^ ratio in our entire sample (~1), indicating ex- 
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treme ISM excitation conditions, with average gas densities of ~(10^-10^) cm"^, and LVG fits 
that keep improving past Tkin=150 K. The 3a upper hmit on r32;$2.54 is certainly compatible 
with such conditions, albeit not imposing any useful additional constraints (high Tgyg because 
of z/sky(3-2) near the 325 GHz atmospheric absorption feature prevented sensitive CO J=3-2 
observations of this system). The high rgs ratio measured for the NW nucleus of this ULIRG 
is actually perfectly compatible with a pure Orion "hot-spot" CO SLED (see 2.4), further 
emphasizing the truly extraordinary levels of ISM excitation. The corresponding molecular 
gas mass is at least Msf(H2)~3.9x10^ Mq (assuming an Eddington-limited SF). A Galactic 
Xco^(3-6)Mq (Kkms~^ p(?)~^ factor is deduced from the LVG solutions with Tkin>65K, 
which corresponds to M(H2)~(4.8-9.6)xlO^M0. Thus MsF/Mtot~0.41-0.81, among the 
highest such fractions for the individualy-studied LIRGs (see Table 1). 



B.IO. Arp55 

This is another clear merger having double nuclei as IRAS 08572+3915, but more widely 
separated (~12", 9.3 kpc) seen in radio continuum maps (Condon et al. 1990) as well as 
optical and CO J=l-0 interferometric maps (Sanders et al. 1988). A tidal tail emerging from 
the eastern nucleus (also the most gas-rich one) is clearly visible in Ha emission (Hattori et al. 
2004). The two nuclei have been separately detected in CO J=3-2 (Leech et al. 2010, Paper 1) 
while the NE nucleus may have been tentatively detected also in CO J=6-5 (Paper 1), with an 
implied high CO(6-5)/(3-2) brightness temperature ratio of r65/32 =0.75 indicating a highly 
excited molecular ISM phase. On the other hand the global HCN/CO J=1~0 brightness 
temperature ratio is ~0.03 (Gao & Solomon 2004), typical of the molecular gas found in 
isolated spiral disks and quite unlike the much higher such ratios in merger/starbursts. 

Using the global r2i and r32, along with the lower limit on the CO/^^CO J=2-l ratio 
(see Paper I) yields a large range of possible average ISM states, all with n~(l-3) xlO^ cm""^ 
and most with surprisingly large Kvir (~70-126). Restricting Kyir to values <40 yields a 
much narrower range of solutions with Tkin=(45-70) K and corresponding Xco~(0.25-0.8) Xi. 
This still considerable Xco range yields Mtot~(2.9-9.2)xlO^M0, while Msf~O.9x1O^M0. 
The latter is ~l/10-l/3 of the total molecular gas mass, leaving ample room for a SF- 
quiescent molecular gas reservoir. In Arp 55 however SF-quiescent is not the cool and mostly 
gravitationally bound gas found in the GMCs of isolated disks. Indeed the large temperatures 
and Kvii- values needed to fit its highly-excited CO lines and R2i^l3 suggest warm and 
strongly unbound gas which, unlike in most other ULIRGs, contains most of the molecular 
gas mass. The low HCN/CO J=l-0 ratio then may not be due to the prevelance of Galactic- 
type SF-quiescent GMCs rather than to a disrupted ISM state that has, momentarily, a 
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low (dense) / (total) molecular gas mass fraction. Short periods during which the dense gas 
reservoir of a LIRG can be strongly disrupted by SF feedback while its IR luminosity remains 
intact are expected in strongly evolving mergers (Loenen 2009). In such a scenario Arp55 
would be one of the few mergers "caught" in a very short act (see also Arp 193). Interestingly 
its HCN(3-2)/(l-0) ratio is <0.3 (Juneau et al. 2009), much lower than in some classical 
(U)LIRGs such as Arp 220 and NGC6240 (Greve et al. 2009). More sensitive ^^CO and 
higher- J HON observations of Arp 55 are necessary to confirm such a tantahzing picture of 
global ISM dynamics. 

B.ll. IRAS 09320+6134 (UGC 05101) 

Ncar-IR NICMOS observations show a single extremely red and unresolved nucleus 
with 6's<0.22" (~170pc), surrounded by a strongly perturbed spiral structure (Scoville at 
el. 2000) with a tidal tail and a large ring seen in the optical (Sanders et al. 1988b). The 
half-light source radius at 2//m is ~450pc while an interferometric CO J=3-2 image reveals 
a gas distribution of ~2"(1.5kpc) in size (Wilson et al. 2008). Its disturbed morphology is 
interpreted as the result of an interaction with another gas-rich spiral (Sanders et al. 1988b) 
and even with a gas-poor dwarf galaxy 17" to the southeast (Majewski et al. 1993). This 
ULIRG contains an AGN luminous in hard X-rays (Imanishi et al. 2003) and is deeply 
dust-enshrouded along most lines of sight, with N(H)>10^'' cm~^ (Imanishi et al. 2001). 

Its CO SLED betrays the presence of a high-excitation gas component already from the 
low-J lines with r2i=1.23 and r32=0.93 (Paper I). These ratios along with the high CO/^^CO 
J=2-l ratio of R2i=18 can be well-fitted by two ranges of LVG solutions, namely: Tkiii=(35- 
65) K, n=3xl03cm-3, K™=22, and Tki„=(100-130) K, n=103cm-^ K^ir=4. The corre- 
sponding average Xco factors are X co~0.4Xi and Xco^O.GXj respectively. We note that a 
small range of solutions with Tkm= (25-30) K, n=10^cm~^ and Kvir=40 also exists but are 
unlikely to be good representation of the average ISM state in this ULIRG since Tkin/Tdust<l 
(considered unhkely, see Paper I) . Moreover their very large Kyir make it even less likely that 
these modest gas temperatures can correspond to such highly unbound dynamical states 
where the molecular gas can be significantly heated by the dissipated supersonic turbulence 
(see Paper 1). Nevertheless adopting this narrow cooler gas solution range would not make 
much difference when it comes to the corresponding Xco factor which is ~0.5Xi. The corre- 
sponding total molecular gas mass then is Mtot~(l-9-2.9)xlO^MQ while Msf~1-6x10^ M©. 

The large fraction of SF gas mass in UGC 05101 (~(55-84)%) is further corroborated by 
the high HCN/CO J=l-0 brightness temperature ratio of ~0.2 (Gao & Solomon 2004) (for 
typical spiral disks this ratio is ~0. 02-0. 03) isofar as HCN J=l-0 is a linear proxy for dense 
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gas mass and most of that gas is involved in star formation. This dense gas reservoir is hkely 
confined into the small (near-IR) -bright nucleus fueling a compact starburst. Such compact, 
~100pc-sized gas-disk/starburst configurations have been revealed for Arp220 (Sakamoto 
et al. 2008; Mathushita et al 2004) where extremely high extinctions (N(H)>10^^ cm"^) 
correspond to significant dust optical depths even at submm wavelengths (Sakamoto et al. 
2010, Papadopoulos et al. 2010a). These can be responsible for Compton-thick AGN, and 
a surpression of high-J CO hues in such systems. This may be the case also for UGC 05101 
as indicated by its low CO(6-5)/(l-0) yet high HCN/CO J=l-0 ratio (see Figure 11, Pa- 
padopoulos et al. 2010b). 

Despite the good one-phase fit obtained for the global CO SLED (and R21) of this 
ULIRG, the observed CO line emission may still be dominated by a phase that does not 
contain the bulk of its molecular gas mass. Indeed, neither the modest densities nor the 
potentially highly unbound average dynamical states implied by the one-phase LVG fits 
would be typical of the HCN-briglit star-forming phase. As discussed in 3.2 the Xco obtained 
from the global low-J CO SLED may thus be suitable only for a gas phase that does not 
contain much of the molecular gas mass. Using the observed Tco/hcn"^ context 
of a 2-phase model along with r^^^^^^-^=2.3A (sec 3.2) yields pco '^•' = 1.136 (Eq. 14). In 
the absence of multi-J HCN observations that could constrain Xhcn we set Xhcn=10Xi, 
the smallest value produced by such studies (e.g. Papadopoulos 2007; Gracia-Carpio et al. 
2008). Thus X^o^ = (rJ5J,/jjc^)-iXHCN-4.3Xi, and for xS=0.5Xi it would be X^o"P''^=2.3Xi, 
which is ~(4-6)x higher than that deduced from one-phase LVG models of the low-J CO 
SLED. Multi-J observations of HCN or other heavy rotor molecules (e.g. CS, HCO+) can 
confirm such potentially high masses of molecular gas at high densities. 



B.12. NGC3310 

The young and intense star formation activity of this UV-bright galaxy is on par with 
M82, and the likely result of a recent merger with a smaller galaxy (Balick & Heckman 
1981; Concehce et al. 2000; Elmegreen et al. 2002). Its gas and dust are warmed and 
disrupted by the presence of the starburst in its central 1-2 kpc. Tidal HI tails, unusually 
large HI velocity dispersions (~40kms~^) for a spiral galaxy (Kregel & Sancisi 2001), and 
large velocity ofi'scts (up to 150kms~^) between molecular cloud and adjacent HII regions 
(Kikumoto et al. 1993) provide further evidence for a dynamically disturbed gas typical of 
mergers. 

This galaxy has been extensively studied by Zhu et al. 2009 as a template of a highly- 
excited ISM, indicated by r2i,32>l ratios and warm IR "colors" (Sloo/^m/S6o^im=l■28). Despite 



-43 - 



its low IR luminosity (L|p^~2.7xl0^'^ L©, Paper I) its global r2i=1.47 and 132=1.21 ratios 
surpass those of Orion cloud SF "hotspots", and are among the highest observed in our 
sample. Its proximity allowed spatially resolved studies of its molecular gas with single-dish 
telescopes obtaining beam-matched observations of the CO(3-2)/(l-0) ratio (see Zhu et al. 
2009 for details). These high global CO line ratios of NGC 3310 can be fitted well by a single 
phase provided Tkm>35K (~Tdust), but remain highly degenerate over the [Tkin, n, Kyir] 
domain up to the highest temperature searched (150 K). Interestingly for Tkin=(35-50) K a 
class of dense (n~(l-3)xl0^cm~^) and strongly unbound (Kvir~20-40) states are possible, 
though only higher- J and/or ^^CO line observations could reduce the wide degeneracy of 
the LVG solutions and designate such phases as those most likely representing the actual 
ISM conditions. Other solution ranges are: Tkin= (70-90) K, n=3xl0^cm~^ and Kvir=22 
while near-virial solutions can be found at high temperatures Tkin=95K, (140-145) K, n=(l- 
3)xl0'^cm^^ and Kvir=4, 2.2. The large degeneracy of possible average ISM conditions 
translates to a significant one for the Xco factors, ranging from ~0.4Xi (for the highly 
unbound states with Kvir=22-40) and up to ~(1.2-2.2)Xi (for the near-virial states with 
Kvir=4, 2.2). 

For Xco=(0.4-2.2) Xi the corresponding Mtot=(0.3-1.5) x 10^ M© while Msf=1. 1 x 10^ cm'^. 
The latter seems to rule out the lowest Xco values while even for the highest ones it is: 
MsF/Mtot~0.73, i.e. most of the molecular gas in NGC 3310 is directly involved in star for- 
mation (for an Eddington- limited SF process) , leaving little room for a massive low-excitation 
molecular gas reservoir. Interestingly the LVG solutions corresponding to the highest Xco 
values are those indicative of a SF phase with n~3x 10^ cm~^, Kvii.=2.2 and large Tkin=(140- 
145) K. We caution though that all this pertains only to the inner ~40" of this galaxy where 
the starburst takes place, a massive molecular gas reservoir with low densities, highly un- 
bound gas motions, and low CO line excitation can still exist beyond it, a situation actually 
encountered in M82 in the form large scale molecular gas outflows (Weiss et al. 2005). 



B.13. IRAS 10565-1-2448 

This is a ULIRG with HII region- type line ratios (Armus, Heckman, & Miley 1989) 
implying an optical spectrum dominated by young stars rather than an AGN. Nevertheless a 
weak contribution from an AGN seems necessary in order to fit its IR/ submm dust continuum 
SED (Farrah et al. 2003). Early r-band (6550 A) imaging has shown this object to be a 
possible triple merger galaxy system (Murphy et al. 1996), while near-IR imaging with the 
HST NICMOS camera shows a luminous primary galaxy interacting with a much fainter one 
~8" (~6.7kpc) to the southeast (Scoville et al. 2000). The primary galaxy is compact with 
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a half-light source diameter at 2//m of ~650pc, and the only source where luminous CO 
J=l-0 line emission is detected (DS98) with a diameter of 1.5" (1.25 kpc). Moreover this 
ULIRG along with VII Zw 031 and Arp 193 arc the only galaxies where kinematic models of 
the CO emission indicate a rotating ring rather than a filled disk gas distribution (DS98). 

Its low-J CO hues (J=l-0, 2-1, 3-2) indicate a well-excited low-J CO SLED (r2i=1.06, 
r32=0.80) expected for a starburst while the large Rio=15 and R2i^l8 are typical for merger 
systems. Two ranges of one-phase LVG solutions that can be found for Tkin/Tdust>l (with 
Tdust=40K) namely: [Tkin, n, K,iJ = [(40 - 50) K, 3 x lO^cm-^,?] and [Tkm, n, K,ir] = [(80 - 
140)K, lO^cm-^A], with the best ones obtained for Tkm=45K and 115 K respectively Both 
also reproduce the CO J=6-5 line luminosity with rg5^''~0.15 (observed value: 0.18±0.055). 

The resulting X^o factors are: Xco=0.75Xi (for Tkin=45K) and 0.60 Xi (for Tkin=115K), 
yielding a total molecular gas mass of Mtot=(3.8-4.80) x 10^ Mq. The minimum such mass 
needed for Eddington-limited SF is Msf=2.6x10^ M©, which is ~54%-68% of Mtot- As in 
UGC 05101 the aforementioned large fraction of molecular gas mass directly involved in 
star formation in IRAS 10565-1-2448 is corroborated by a high HCN/CO J=l-0 brightness 
temperature ratio of rHCN/co~0.17 (Gao & Solomon 2004). However all the 1-phase LVG 
solutions obtained are hardly representative of a dominant dense and self-gravitating gas 
phase but more typical of the "cloud-envelope" warm and diffuse gas often found in the ISM 
of mergers. Moreover all one-phase LVG fits fail to reproduce the R2i^l8 lower limit, with 
R-2i'^''=12 being the highest vahic obtained (for the Tki„=115K solution). A 2-phase model 
with an assumed (h)-phase alcviatcs this and yields an (l)-phasc with Tkin=65 K, n=10^ cm~^, 
and Kvir=13 that has R2i~20 and Xco=0.4Xi. The computed Xco ^''■'~0.75Xi is within the 
range of the one-phase values, a result of the low xS, the large luminosity contribution of 
the (l)-phase (pco~'^^~4.5), and the modest X^o^=2.2Xi of the (h)-phase (see 2.4). 

Setting constraints on the (h)-phase and Xco'' by using the observed CO J=3-2, 6-5 
and ^^CO J=l-0 lines (assumed emanating mostly from that phase), does not change the 
aforementioned picture, at least when it comes to the values of Xco "^^^ ■ Indeed while the 
corresponding LVG models do recover dense gas solutions (n~10^cm~^), substantial Kyi^ 
values (~13) keep the corresponding X^o^ well below Galactic. In this model the (l)-phase 
contribution, and thus the influence of xS on the xi^P'^\ is negligible. We must stress 
however that other high-J CO line observations, and crucially multi-J HCN observations 
are necessary to confirm this picture. If luminous such lines were to be discovered for this 
ULIRG, they could substantially raise Xco^ and thus Xco"^'^-' as discussed in 3.2. 
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B.14. Arp299 

This is a spectacular merger of two galaxies, IC694 and NGC3690, with luminous CO 
J=l-0 emission detected in the nuclei of both (which are 22" - 4.4 kpc, apart), and the 
interface between them (Sargent et al. 1987; Sargent & Scoville 1991). A long HI tidal tail is 
also seen extending out to ~180kpc (Hibbard & Yun 1999). Along with the Antennae galaxy 
this (U)LIRG has been considered an early template for mergers, documenting their ability 
to rapidly funnel the gas in the disks of the progenitors into compact regions (Sargent & 
Scoville 1991). Indeed the two nuclei in Arp299 contain ~80% of its molecular gas reservoir, 
fueling intense star-formation, while the nucleus of IC 694 contains also an AGN (Sargent & 
Scoville 1991 and references therein), and is the primary source of the bolometric luminosity 
of this spectacular system (Charmandaris et al. 2002). The most detailed molecular line 
study of this template merger system is the interferometric study of CO, ^^CO and HCN 
J=l-0 hne emission by Aalto et al. 1997 which finds most of the HCN-bright emission 
located in compact regions (<310pc) in the two nuclei of the merging galaxies, and unusually 
large CO/^^CO line ratio variations (from ~60 in the IC694 nucleus, to ~5-10 in its disk). 
Actually the state of the gas in the IC 694 nucleus with the largest HCN/CO and CO/^^CO 
ratios exemplifies the effects of highly turbulent environments found in mergers, with much 
of the molecular gas mass "resettled" at high densities (boosting the HCN/CO ratio) while 
the violent disruption of GMCs creates a diffuse warm and highly unbound "envelope" phase 
dominating the low- J CO lines (boosting the CO/^^CO ratios). 

While we have no CO J=6-5 measurements for this galaxy (too extended for the narrow 
JCMT beam at 690 GHz) our CO J=3-2 and J=4-3 measurements indicate high CO line 
excitation. A one-phase LVG model of the global CO (3-2)/(l-0), (4-3)/(l-0) and CO/^^CO 
J=l-0 line ratios gives the best fit for Tkin=30K, n=10^cm~^ and Kvir=40 which indicates 
strongly unbound gas motions. Other average gas states are also possible (e.g. Tkin=40 K, 
n=3xl0^cm~^ and Kvir=22) but all have Kvir^22. The corresponding Xco factor over the 
LVG solution range is Xco=(0. 35-0.42) Xi, with the largest value corresponding to the dense 
gas solution. These yield Mtot~(l-l-2) x 10^ while Msf~1.9x 10^ M0, which clearly favors 
the higher Xco values and indicates that ~100% of the molecular gas mass is associated with 
SF sites. Larger values of Xco (and Mtot) remain possible however if most of the molecular 
gas mass in Arp 299 is contained in the HCN-bright phase with large Xhcn factors. Multi-J 
CO, ^^CO and HCN, H^^CN line imaging would be ideal in determining this as a function of 
position within a template merger whose ISM conditions seem to encompass the full range 
possible, from quiescent disks to compact starburst nuclei. 



B.15. IRAS 12112+0305 



This ULIRG consists of a strongly interacting pair of galaxies whose nuclei are ~2.9" 
(~4kpc) apart (Carico et al. 1990; Scovile et al. 2000), while tidal tails and "plumes" are 
visible both in optical (Surace et al. 2000) and near-IR (Scoville et al. 2000) wavelenghts. 
CO J=l-0 interferometry finds the bulk (~75%) of the emission emanating from the NE 
nucleus of this system (Evans et al. 2002), which is unresolved with a size of <2" (<2.7kpc), 
and also the only one tentatively detected in CO J=6-5 (see Paper I) . Thus our one- phase 
LVG radiative transfer modeling pertains only to the NE nucleus of this strongly interacting 
system. The observed r32=1.58 ratio is the second highest in our sample (the highest one 
found in another double nuclei system: IRAS 22491-1808) and indicates the presence of high- 
excitation molecular gas. This surpasses the corresponding ratio of the SF "hot-spots" in the 
Orion molecular cloud, while the r32>r2i inequality is another indication of a qualitatively 
different CO SLED, irreducible to a mixture of SF and SF-quiescent gas (see Paper I). 

All one-phase LVG solutions for IRAS 12112+0305 indicate average gas states with 
n~(3x 10^-10^) cm~^. Nevertheless lack of ^^CO line observations leaves serious degeneracies 
with e.g. [Tkin(K),n(cm-3),Kvir] = [60, 10^ 13], [125, 10^ 1] being nearly equivalent. The 
corresponding Xco~(1.5-5.3) Xi, with Galactic values of ^^(4-5) Xi obtained for the warmest 
(—(120-150) K) and densest (=(l(P-l(fi) cwT^) states. The presense of such high densities 
is corroborated by a high HCN/CO J=l-0 ratio of rHCN/co=0.16 (HCN from Gracia-Carpio 
et al. 2008, CO from Paper I), though without offering any constraints on the temperatures 
of the dense gas or its all-important kinematic state. However even the warmest and densest 
one-phase solutions remain rather poor fits of the "hot" r32~1.6 and comparatively "cooler" 
r2i~0.9 ratio. Indicatively such a solution (n=10^cm~^) with Tkin=150K and Kvir=4 yields 
r32~r2i~l as the high density, and modest K^ir act to keep the CO lines well-thermalized 
but also optically thick up to high-J levels. 

A two-phase model of the global CO line ratios in this intriguing ULIRG cannot use 
an Orion-derived (h)-phase CO SLED as the observed r32>r32'' (Orion). Using only the 
CO (3-2)/(2-l) ratio r32/2i = 1.76 ± 0.48 to constrain the high-excitation (h)-phase yields: 
n=3xl0^cm-3, K™=22, and Tkin>140K (for Tkin=140K, r^2/2i=l-28-r32/2i-(T), with a fit 
that keeps improving past Tkin=200K (where r32/2i=1.41). This is expected since r32/2i>l 
indicates a thermalized J=3-2 line (and thus n>nCTit,3-2~10^ cm~^), with low/modest opti- 
cal depths (ensured in part by the large Kyir), anc? warm enough gas to populate the J=3 level. 
The molecular gas mass of this phase can be estimated from LQQ^^_Q^=(l/r32'*)LQQ^3_2)~4xl0^ Li 

(where r32^~4) and the computed Xco^~ 1 Xi, thus Mh_ex(H2)-^4x 10^ M0, and similar to that 
needed to fuel an Eddington-limited SF in this ULIRG of Msf(H2)~5x10'^ M©. For the 
(l)-phase L^Q^_Q^~6.1xlO^ Li luminosity there are no constraints on its corresponding Xco 
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factor. Setting X^o^=(0.5-5) Xi yields X^o ^ ^=(0.7-3.4) Xi and a corresponding Mtot~(7.1- 
34.3)xlO^M0. Multi-J HCN and R'^^CN as well as ^^CO line observations are critical for 
reducing the degeneracies of the radiative transfer models for the dense and the low-excitation 
gas, and their Xco factors. 



B.16. Mrk231 

This is an archetypal ULIRG/QSO and the most IR luminous galaxy in the Revised 
IRAS Bright Galaxy Survey (RBGS; Sanders et al. 2003). It has a compact nucleus, sur- 
rounded by irregular incomplete rings of star formation along with a small tidal arm that 
contains numerous blue star-forming 'knots' (Farrah ct al. 2003 and references therein). It 
has large amounts of molecular gas in a compact (0.85", 700 pc) nearly face-on disk (Bryant 
& Scoville 1996; DS98) which allows nearly unobscured view towards an optically lumi- 
nous AGN which classifies this object also as a Seyfert 1 galaxy. Mrk231 was also the first 
ULIRG for which high-J CO fines (J=4-3, 6-5) were detected (Papadopoulos et al. 2007) 
while recent SPIRE/FTS observations with the HSO revealed luminous high-J CO fines up to 
J=13-12 emanating from AGN-induced X-ray Dissociated Regions (XDRs) (van der Werf et 
al. 2010). Its well-excited high-J CO fines, the largest HCN/CO J=l-0 ratio (rHCN/co~0.29) 
observed among ULIRGs (Papadopoulos et al. 2007), and a very high CO/^^CO J=2-l ratio 
exemplify the extraordinary average molecular gas states possible in such systems (similar 
to the IC694 core in the Arp299 merger system). Finally the availability of HCN J=4-3, 
1-0 measurements along with the high-J CO lines allowed a detailed study of its dense gas 
phase (Papadopoulos et al. 2007), and make Mrk231 a good case study for the effects of the 
large M(n>10^cm~"^)/Mtot fractions on the global Xco of ULIRGs. 

A one-phase LVG fit predictably fails to converge on any average state as it cannot 
accomodate both the high densities needed to excite the CO J=3-2, 4-3, and 6-5 fines and 
maintain the very low optical depths needed for reproducing the very large Rio=47 ratio. 
Indicatively the best one-phase LVG solutions: Tkin=(115-150) K, n~(300-10^) cm~^ and 
Kvir=7-40 yield r^\''^^~0.83-l.l (obs=0.89±0.12), rj2''^^~0.5-0.7 (obs=0.72±0.13), ri5'^^~0.22- 
0.3 (obs=0.8±0.2) and r^5^^~0.03 (obs=0.42±0.12, but not used for the LVG fit), while 
R^^g) =25-32 (obs=47±16). Thus the average ISM state deduced by a one-phase LVG model 
can be responsible for up to J=3-2 (and the ^^CO J=2-l) line emission but severely under- 
prcdicts the CO J=4-3, and 6-5 line luminosities. The corresponding Xco~(0.25-0.45) Xi is 
low but within what was thought as the ULlRG-appropriate range (for Mrk231 DS98 give 
Xeo~(0.7-0.8)Xi). These yield Mtot~(1.8-3.2)xlO9M0, but a minimum Msf=5.5x109 
clearly favors the higher Mtot (and Xco) among those computed with the one-phase model. 
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For a 2-phase model we use the results by Papadopoulos et al. 2007 of warm (~(85- 
140) K), diffuse gas (~300cm~^) that is highly unbound (Kvir^l), and dense gas (~(1- 
3)xl0^cm-3, Tkin=(40-70) K) responsible for the bright CO J=6-5, HCN J=l-0, 4-3 lines. 
The latter phase is virialized with Xhcn~(10-25) Xi (for our numerical factor in Equation 
3 rather than the one in Papadopoulos et al. 2007) and t^^q^^^j^^2-2.54 (high value cor- 
responding to the lower Xhcn values). Thus Xco^ = (l/r[!)^yjj(^j^)XHCN~(3.94-12.5) Xi while 
for the diffuse and HCN-dark phase: xS~(0.4-l) Xj. For an Tqq^^^j^—3.4:5 we then find 
pi'o"^^~0.36-0.73 (yielding xio^ = (3.94-12.5) Xj), and an X^o"''^^~(3-7) Xi (for the smallest 
Xco). Thus a Galactic Xco factor applies in Mrk231, a result of a large mass of dense gas 
with high xio'' values, themselves a result of the high gas densities and their virial dynami- 
cal states. 



B.17. Arpl93 

This galaxy is a template of a young merger-induced starburst with an age of fewx 10^ yr 
and near-lR imaging showing a highly dust-enshrouded nucleus, and tidal tails (Smith et 
al. 1995). A disturbed highly inclined disk is also seen in NICMOS near-IR imaging with 
reddening increasing towards the NW within the disk (Scoville et al. 2000), and a similar 
extinction distribution implied by a CO J=l-0 interferometer map (Bryant & Scoville 1999). 
Extended radio continuum emission (~3"x3.7") is seen at cm wavelengths with two nuclei 
~1" apart and a spectral index of q;~0.6 {Si,(xu~°'), consistent with SF- related non-thermal 
synchrotron emission (Clements & Alexander 2004). The latter study also finds significant 
decoupling of the HI and CO-bright H2 gas motions and distributions. The young age of 
this disk-disk merger/starburst is corroborated by several extra-planar highly luminous star 
clusters (Scoville et al. 2000) and many bright HII regions (Clements & Alexander 2004). 

High S/N CO J=l-0 and 2-1 interferometric imaging confirmed a highly inclined 
2.8"x0.8" (~1.3kpcx0.37kpc) gas disk whose CO emission is best fit by a molecular ring 
(DS98). A compact southeast CO source coinciding with the secondary cm continuum emis- 
sion peak in Arp 193 has the properties of a giant molecular core, and is much warmer and 
denser than the rest of the disk. Along with a region in Mrk273, and the western nucleus of 
Arp 220 (DS98), it is one of the most extreme SF regions in the local Universe, forming stars 
very near or at the Eddington-limit set by the radiation pressure from young massive stars 
on the dust of the dense molecular gas (i.e. LiR/Mdense(SF-region)~(250-500) L0/Mq). One- 
phase LVG solutions to its CO line ratios yield low-density solutions with n~(l-3)xl0^cm~^ 
but large temperature degeneracies with Tkm=40 K, (60-75) K, (80-115) K and corresponding 
Kvir=22, 4 and 13 values. Warmer solutions with Tkin=(120-150) K and even larger Kvir(=40) 
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values also exist. Considering only the gas states with Kvir<20 yields Xco~(0.35-0.75) Xi and 
Mtot~(l-6-3.5) xlO^ Mq, with the highest X^o (and Mtot) values corresponding to the states 
with the smallest K^i^- The minimum molecular gas mass for an Eddington- limited star 
formation is MgF~0.9xl09 Mq. 

A young starburst and the strong SF feedback expected in such extreme SF events 
could be responsible for Arp 193 being the only LIRG where, despite a large SF-powered 
IR luminosity (~2.2xlO"L0) and HCN/CO J=l-0 ratio (~0.044, with HCN J=l-0 from 
Carpio-et al. 2008, and CO J=l-0 from Paper I), the average state of its molecular ISM, 
is compatible with a complete absence of a dense molecular gas reservoir (Papadopoulos 
2007). This is indicated by its low HCN (4-3)/(l-0) and (3-2)/(l-0) line ratios of r43<0.12 
(3(t) and r32=0.22±0.03 (Papadopoulos 2007; Gracia-Carpio et al. 2008). Using them as 
constraints yields a narrow temperature range of dense gas solutions with Tkin=(50-60) K, 
n=3x 10"^ cm~'^ and Kvir=14 (rather large for a dense gas phase). The widest range of solu- 
tions is at: Tkin=(65-95) K, n=3xl0^cm^^ and Kvir~l. Unlike the former set of solutions 
the latter reproduces also the CO(3-2)/HCN(l-0) line ratio of ~17±4 (which can be reason- 



ably assumed as an additional constraint on the denser SF gas phasci). Using the upper 



limit of the HCN(4-3)/(l-0) ratio as a detection (as to obtain the densest possible LVG so- 
lutions) only pushes the aforementioned set of dense gas solutions to Tkin= (90-105) K, while 
the lower density /virial ones are now found at Tkin=(110-150) K (both sets of solutions can 
now reproduce the CO(3-2)/HCN(l-0) ratio). 

For a 2-phase model where the HCN-bright gas is set as the (h)-phase we obtain 
Xco '''^''~(0.70-1.75) Xi (with the low values corresponding to the dense/unbound solutions). 
This is is ~(0.9-2)xXco(DS98), and yields Mtot~(3.3-8.1)xlO9M0. Thus Arp 193 stands as 
a case where multi-J line observations of high-density tracers such as HCN do not signifi- 
cantly change the picture obtained by the low-J CO lines regarding the global Xco factor. 
This may be because this LIRG is a rare merger, "caught" in the act of having momentarily 
dispersed/consumed its dense gas mass while still having large SF-powered IR luminosities, 
a brief state of affairs expected because of SF-feedback during the evolution of merger star- 
bursts (Loenen 2009). We note however high- resolution CO imaging revealing the "excess" 
CO J=l-0 line emission of the (l)-phase to be a cold SF-quiescent phase can still push the 
global Xco towards near-Galactic values via a higher xS. 




^thc non-dctcction of CO J=6-5 in this LIRG is likely due to a pointing offset of the JCMT 
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B.18. NGC5135 

Optical and near-IR i/ST observations reveal this Sy2 galaxy to be a grand design spiral 
with star formation across its arms but being most intense in its central ~lkpc (Martini et 
al. 2003). The inner 1.6 kpc is also where most of its cm radio continuum emerges (Condon 
et al. 1996). Optical and UV studies by Gonzalez Delgado et al. (1998) reveal a young 
nuclear star burst with gas hkely provided by a bar instability. Its HCN/CO J=l-0 hne ratio 
of ~0.087 (Gao & Solomon 2004), while not as high as in ULIRGs, it is nevertheless ~(3- 
4)x higher than in spiral disks. This indicates a significant dense gas mass fraction, even 
in the absence of a strong merger and/or galaxy interaction. Lack of the latter makes the 
large Rio=26 and Rio=13 measured in NGC5135 an exception as these are typical mostly in 
mergers. Interestingly this galaxy is known for powerful gas outflows of hundreds of km/s, 
driven by its young starburst (Gonzalez Delgado et al. 1988). 

One-phase models of its CO line ratios yield [Tkin, n, K^ir] = [(30-40)K, 3 x 10^ cm~^ 22], 
with warmer (but also more degenerate) solutions such as: [(60 — 110)K, 10^ cm^'^, 13] also re- 
maining possible. For the entire range of solutions the corresponding Xco~(0.35-0.45) Xi, sim- 
ilar to Xco"'^'' (LTE) for Tkin=(30-70) K (the temperatures of the best solutions), as expected 
for the low/modest optical depths of the LVG solutions. For these Xco values Mtot=(l-l- 
1.4)xlO^Mo while Msf^SxIO^Mq, which corresponds to ~20%-30% of the total molecular 
gas mass. Nevertheless, with the distribution of the molecular gas unknown, a Galactic Xco 
for disk-distributed gas not intimately involved with star-formation remains a distinct pos- 
sibility. In such a case the true Xco and Mtot can easily be ~(10-15) times higher, reducing 
Msp/Mtot to only few%, and more typical of spiral disks. As described in 3.4 CO line imaging 
observations rather than global CO SLEDs can be used to determine the existence of such 
a cold molecular gas disk. 



B.19. Mrk273 

This extraordinary ULIRG shows up as a double-nuclei system in NICMOS near-IR 
observations where the nuclei are only ~1" (740 pc) apart (Scoville et al. 2000). A strikingly 
long tidal tail (~1', 44 kpc) as well as a decoupling of stellar and gas motions (Tacconi 
et al. 2002) are the telltale signatures of a strong merger. Its northern nucleus contains 
a compact hard X-ray luminous Sy2 AGN (Xia et al. 2002), but the dust emission SED 
of this ULIRG is dominated by a starburst (Farrah et al. 2003). This galaxy actually 
contains the most extreme starburst region among all ULIRGs imaged in the DS98 study, 
with an IR luminosity of ~6.7xlO^^L0 emanating from its core of only ~0.35"x(<0.2") 
(260pcx(<149pc)) diameter. This makes the compact nuclear CO hne source of Mrk273 
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akin to the west nucleus of Arp 220, and along with Arp 193 these three ULIRGs host the most 
prodigious SF events in the local Universe (DS98). For the core of Mrk273 the IR brightness 
is: (TiR~LiR/(7rR^Q)~10^'^ kpc^^, i.e. the Eddington limit for radiation-pressure-regulated 
star formation (Thompson 2009). In the framework of an Eddington-limited SF this makes 
the compact nuclear starburst in Mrk273 a maximal SF event, and its astounding CO 
linewidth of ~1060kms"^ (DS98) may be a direct outcome of radiation pressure becoming 
dynamically important as to affect bulk gas motions. 

The global CO line ratios, while typical of a star-forming ISM they are not typical of 
an extreme starburst as Mrk273 certainly is. The R2i=7 ± 2 ratio in particular is low, 
found mostly for GMCs in SF-quiescent spirals (R2i>20 is observed for extreme starbursts). 
The compact starburst core of Mrk 273 is surrounded by extended CO hne emission (~3"- 
7", i.e. (2.2-5.1) kpc, DS98) where ordinary star- forming and even SF-quiescent ISM could 
be containing significant fractions of the total molecular gas mass. The optimum LVG 
solutions obtained for the global CO ratios are: [Tkm, n, Kvir] = [30-70,3x 10^, 1]. Apart from 
the elevated temperatures such conditions are typical of Galactic GMCs, while the Xco~2 Xi 
of the best LVG solution (at Tkin=45K) is higher than the so-called ULIRG value (~0.8Xi). 
The corresponding Mtot~lO^°M0 while a minimum mass of Msf~3.3x10^ M© is needed to 
fuel an Eddington-limited star formation in this system. 

The total molecular gas mass of Mrk 273 can be higher still if the mass contributions 
of a colder low-dcnsity and SF-quiescent phase, and a warmer denser SF one are accounted 
separately. A clear indication of a highly-excited SF gas phase is given by the CO J=6-5 
line which, while not as bright as in some other (U)LIRGs (e.g. Mrk231, NGC6240), yields 
a r65 that is ~5x higher than that anticipated from the best LVG fit of the lower- J CO fines. 
A high HCN/CO J=l-0 ratio of rHCN/co=0.14 (using the CO J=l-0 value from Paper I 
and the HCN J=l-0 from Carpio et al. 2008) and a well-excited HCN J=3-2 fine with 
r32(HCN)=0.49±0.11 (Gracia-Carpio et al. 2008) offer more evidence for the presence of a 
dense gas phase. A two-phase decomposition using only the CO 3-2, 6-5 and ^^CO J=2-l 
lines to constrain the high-excitation phase yields Mh-ex(Il2)~6x 10^ M0, and an effective 
xio"''''^~2.5Xi. The latter yields Mtot(H2) = 1.3xlO^°M0, of which ~53% correspond to the 
low CO brightness (l)-phase. Thus despite its spectacular starburst, Mrk 273 is an example 
of a ULIRG that contains also large amounts of cooler low CO-excitation gas. Use of the 
two available HCN lines to determine the propreties of the (h)-phase does not change this 
picture, and can actually raise Xco~^^^ up to ~(4-5.5) Xi. 
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B.20. 3C293 

This powerful radio galaxy was the first F-RII source (i.e with edge-brightened radio 
lobes, see Fanaroff & Riley 1974) to be detected in CO locally (z=0.046), revealing large 
amounts of molecular gas distributed in a rotating disk of 1" (~6.2kpc) centered in its nu- 
cleus, and with a large CO hnewidth of ~900kms~^ (Evans et al. 1999). Its disturbed 
morphology maybe due to an interaction with a companion galaxy ~40"(35kpc) to the 
southwest, while a powerful optical/IR-luminous jet interacts with the ambient ISM (Floyd 
et al. 2006). The SF efficiency of the galaxy is low (^SLq/M©) and along with a SFR~(6- 
7)M0yr^^ it is typical of ordinary spirals. The modest global CO r2i=0.74 and r32= 0.44 
ratios (Paper I) are also consistent with a SF-quiescent ISM. Thus the high excitation levels 
of the CO J=4-3, and 6-5 discovered in this system were certainly a surprise. Their inter- 
pretation as due to an AGN-driven jet-ISM interaction provided the first known example 
where "mechanical" AGN feedback globally affects galaxy-size molecular gas reservoirs (Pa- 
padopoulos et al. 2008, 2010b), though similar examples with lower-power AGN and much 
smaller molecular gas reservoirs have been reported much earher (Matsushita et al. 2004). 

The physical conditions compatible with the luminous CO J=4-3, 6-5 lines were studied 
in detail by Papadopoulos et al. 2010b (see their Table 4) where dense (~3x (10^-10^) cm~^), 
warm Tkin~ (75-300) K and strongly unbound Kvir~ 15-50 states have been deduced. The 
corresponding X^o^~(0.8-1.5) Xi, with computed r£^=2-3.2 and LS5^'=(l/r£^)xL43=(1.46- 
2.33)xlO^Li. This AGN-excited gas will have Mh_cx~(l-l-3.5) xlO'^ M©, much higher than 
the MsF~5xlO''M0 of warm and dense gas expected fueling its star formation. Setting a 
Galactic Xco=5Xi for the remaining (l)-phase CO J=l-0 luminosity of the SF-quiescent disk 
phase yields an effective xio~'''^'' = (2.9-3.9) Xi, which are near-Galactic values. Indicatively, 
a one-phase fit of only the CO 1-0, 2-1, and 3-2 lines (any LVG fit including the CO J=4-3 
and/or J=6-5 predictably fails) yields maximum values of Xco=(1.2-1.6) Xi, obtained for the 
solutions with the lowest Kyir values (~2-4), but values as low as Xco~0.5Xi also probable. 



B.21. IRAS 14348-1447 

This is an impressive ULIRG/merger of two gas-rich spirals with their two nuclei 3.4" 
(5.2 kpc) apart, and clear indications of tidal features and strongly disrupted disks seen 
in EST optical and near-IR NICMOS images (Evans et al. 2000; Scoville et al. 2000). 
Interferometric imaging of CO J=l-0 has shown that both nuclei contain copious amounts 
of molecular gas but remain unresolved with size <2.5" (3.8 kpc) (Evans et al. 2000). Using 
radio continuum measurements to obtain size estimates for the two star-forming cores of 
this ULIRG yields very compact regions with diameters of D~200pc (Condon et al. 1991), 
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similar to the compact molecular gas disks found in Arp220 and Mrk231 (DS98). 

The observed CO J=l-0, 2-1, 3-2 and ^=^C0 J=2-l lines (Paper I) are compatible with 
a wide range of conditions having Tkin>35K with the best ones found for Tkin=(40-85) K, 
n=10^cm~^, and Kvir=4. Quite unlike most merger/ULIRGs, all good LVG solutions corre- 
spond to Kvir~l-4. Nevertheless Xco~(0.65-1) Xi, i.e. much lower than Galactic. The corre- 
sponding Mtot~(l.l-1.74)xlOi°M0 while the minimum mass needed for Eddington-limited 
star-formation is Msf~5.6x10^ Mq, which amounts to ~30%-50% of the total molecular 
gas mass. In IRAS 14348-1447, as in other ULIRGs, such high mass fractions of what is 
presumably a dense and warm gas phase arc not reflected by the global CO SLED and 
^^CO lines available, dominated by lower-density gas. Unfortunately neither high-J CO nor 
HCN line observations exist for this ULIRG to set independent constraints on the mass 
and physical properties of its dense SF molecular gas. A 2-phase model with an assumed 
(h)-phase CO SLED (section 2.4), while not necessary in this case (as a wide range of good 
one-phase fits exist), can be used to assess the impact of a massive dense gas phase on Xco- 
Apart from n<10'^cm~^ the (l)-phase remains essentially unconstrained (as expected), with 
xS-(0.5-1.12)Xi. For the computed p^o"^^=0.85 this yields X^o"''''^ = (1.4-1.7) Xj. 



B.22. Zw 049.057 

This galaxy has a SF-powered l{^^~1.2 x 10^^ Lq typical for SF spiral disks (e.g. NGC 7469) 
and is one of the lowest IR luminosity galaxies in the IRAS BGS. NICMOS near-IR imaging 
with the HST indeed show a highly inclined disk-dominated system with a heavily dust- 
enshrouded nucleus (Scoville et al. 2000). Despite its disk-like appearance its near-IR 
emission is better flttcd by a r^^^ rather than an exponential disk proflle (Scoville et al. 
2000), while its Ljr/Lb ratio is almost as extreme as the much more IR-luminous Arp220 
(Planesas et al. 1991). Early interferometric images of CO J=l-0 found large molecular gas 
mass within a 1.3kpc-sized region (Planesas et al. 1991, for the cosmology adopted here). 

Its CO SLED remains well excited up to J=4-3 and also has substantial J=6-5 emission 
(see Paper I) while both its ^^CO J=l-0, 2-1 lines are detected and yield Rio~16 and R2i~24, 
typical of merger-driven starbursts rather than disk-dominated LIRGs. One-phase LVG 
models reproduce all CO, ^^CO hue ratios up to CO J=4-3 though they still underpredict the 
luminosity of the latter by ~ 1.2(7. Unfortunately, even with this large set of lines significant 
degeneracies remain. Indicatively the two best solutions: [Tkm, n, Kvir]=[35, 3 x 10^,22], 
[105, lO'l 13] yield CO line ratios: [rsi, rga, r43, rgs, Rio, R2i]ivg=[l-05, 0.75, 0.4L 0.05, 22, 21], 
[1.05, 0.76, 0.40, 0.02, 21, 16]. Much warmer gas with [Tkm, n, Kvir] = [130 - 140,10^4] is 
also possible and yields: [r2i, r32, r43, rgs, Rio, R2i]ive=[l-06, 0.86, 0.60, 0.17, 14, 13], now in 
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better agreement with the observed res (even if not used in the fit) but a rather worse fit of 
R21. Nevertheless the corresponding X^o values span a rather narrow range ~(0.35-0.6)Xi, 
yielding Mtot~(3-5.2)xlO^M0. This is comparable to the minimum molecular gas mass 
required for Eddington-limited star formation of Msf(H2)~4.7x 10^ M©, yet all the densities 
deduced from the LVG solutions of the global CO SLED are far from being representative 
of what is typically a much denser SF phase with small/modest Kvir values. 

A two-phase decomposition of the average ISM state using the CO J=3-2, 4-3, 6-5 and 
the ^•^CO transitions to determine the (h)-state converges to [Tkin, n, Kvii.] = [110-140,10'^,4], 
with Tkm=140K corresponding to the best solution (and Xco''=0.6 Xi). Higher density so- 
lutions with n=10"^cm~^ exist over Tkin=(40-65) K but represent poorer fits and correspond 
to unbound states with Kvir~13. This is a general aspect of all (h)-phase LVG solutions 
(e.g. even if only the CO J=4-3, 6-5 and ^^CO lines are used), namely the best fits occur 
at higher temperatures, lower densities, and more gravitationaly bound states. High density 
solutions (~10^cm~^) and thus more typical for the presumably star forming (h)-phase are 
found at lower temperatures (~(35-70)K), and strongly unbound states (Kvir~13-40), but 
yield poorer fits. Moreover the high Kyij. of the denser/cooler LVG solutions yield similarly 
low X£^ (e.g. for Tkin=50K, n=10^cm-3 and Kvir-13 it is x£^=0.7Xi). 

It is the combination of weak ^^CO J=l-0, 2-1 and well-excited high-J CO lines that 
nevertheless still have rj+ij<l which forces LVG solutions to a parameter space with low 
CO line optical depths and modest gas densities. This is why K^ij. is high towards the 
low-Tkin/high-n domain (as to keep CO line optical depths low), and why the need for high 
Kvir's is relaxed towards higher Tkin which now becomes partly responsible for the low optical 
depths. Average gas densities towards the high Tkin regime cannot be high as this would 
yield rj+ij>l in the optically thin domain. If we omit the CO J=6-5 line from the fit (as its 
luminosity can be highly uncertain) and also leave out the ^^CO J=l-0 line (which may have 
a larger contribution from a cooler disk phase) we obtain LVG solutions with densities of 
~(l-3)xl0^cm~^ that are more typical of SF gas, yet with high enough Kvir~13-22 that still 
yield a low corresponding Xco''~(0.65-0.75) Xi. Moreover, with an essentially unconstrained 
(l)-phase, and no evidence for a cold disk, a 2-phase gas would yield Xco '''^'* similar to that 
obtained from 1-phase model. 

Thus the molecular gas reservoir in this LIRG is expected to be dominated by its SF 
component, yet neither the average densities nor the strongly unbound dynamical states 
deduced from its available global CO SLED are typical of the latter. If this is confirmed by 
future HCN or other heavy rotor molecular line observations (or more high-J CO and ^^CO 
transitions) it may come to resemble the average ISM state deduced for Arp 193. 
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B.23. Arp 220 and NGC 6240 

These are the two (U)LIRGs in our sample that currently have the best molecular line 
data tracing their dense gas (i.e. muIti-J HCN, CS, HCO+ lines), which are used in a 
detailed study of its physical conditions (Greve et al. 2009). The presence of a massive 
dense gas phase in both of them is well-established with large HCN/CO J=l-0 line ratios 
(~0.19(Arp220), ~0.08(NGC 6240)), well-excited HCN and CS higher-J lines, yet also the 
large CO/^^CO hne ratios typical for merger-driven extreme starbursts (~43-45 for J=l-0). 
The latter implies very low CO line optical depths, and thus runs counter what is expected 
for dominant high-density molecular gas reservoir (rj+i j(CO)^l), necessitating the 2-phase 
models used in such systems (Aalto et al. 1995; Papadopoulos & Seaquist 1998; this work). 

Predictably one-phase models of the CO and ^^CO hnes (see Paper I) fail to converge on 
any good set of solutions for both galaxies, yielding [Tkm, n, Kvir] = [(30 — 90) K, 10'^ cm~^, 40] 
(Arp 220) and [Tkm, n, K^ir] = [(15 - 40) K, 3 x 103cm-3,40] (NGC6240). For Arp 220 even 
the CO J=l-0, 2-1, 3-2 lines are impossible to fit with one phase since r2i=0.67±0.07 (see 
Paper I) indicates a diffuse gas phase that leaves even the J=2-l line subthermally excited, 
while r2i=0.97±0.14 marks the emergence of another, denser and warmer phase from J=3-2 
and higher. The corresponding one-phase Xco~0.3Xi for both systems, a low value due to 
the very large Kyir values "forced" by the extreme large CO/^^CO hne ratios obtained for 
these systems. This is comparable to X^o*^'"^ (LTE) (for the Tkin=(30-40) K range of one-phase 
LVG solutions), as expected for optically thin CO hne emission. 

This low Xco factor yields Mtot~1.85xl09 M0(Arp22O) and ^2.5x10^ Mq (NGC 6240), 
while the minimum molecular gas mass needed for fueling Eddington-limited SF in these ex- 
treme merger/starbursts is Msf~4x10'' MQ(Arp220) and ^1.5x10^^ Mq(NGC 6240). Given 
that MsF/Mtot^O.6, and in the case of Arp 220 actually >1, it becomes obvious, even in the 
absence of heavy rotor molecular line observations that the true Xco factor for these two galax- 
ies must be higher. From the analysis detailed in 3.2.1 the multi-J HCN LVG solutions in a 
2-phase LVG model yield: xio~P^'^~(2.4-4.5) Xi(Arp 220) and xio"''''^~(l-3.3) Xi(NGC 6240) 
with the high values corresponding to virial dynamical states for the dense gas. The new 
Mtot estimates are: ~(1.5-2.8)xlOi°M0(Arp22O) and ~(O.8-2.8)xlO^°M0(NGC624O). In 
the case of Arp 220 these amount to ~40%-77% of the dynamical mass of its molecular gas 
disks (Mdyn from DS98), while for NGC 6240 it surpasses the total mass contained within its 
central radius of r<l" (480 pc) and argues for significant molecular gas mass corresponding 
to the CO line emission seen beyond its central CO J=2-l peak (see Tacconi et al. 1999 
for CO J=2-l for a high resolution interferometry map). In the latter case the deduced 
molecular gas mass may actually dominate the mass distribution out to much larger radii 
than deduced by Tacconi et al. 1999. Multi-J HCN line interferometric maps, along with at 



-56- 



least one H^^CN transition, at resolutions of <1", will be valuable in determining the Xhcn, 
the effective and X^o~^^^ of these template merger systems. 

B.24. IRAS 17208-0014 

This classic ULIRG is the most luminous in our sample, it has a compact gas-rich 
nucleus (Planesas et al. 1991; DS98), a very disrupted disk, and tidal tails indicating a 
merger (Murphy et al. 1996; Scoville et al. 2000). In the interferometric study by DS98 this 
is the ULIRG whose compact gas disk (~1.4kpc) has the largest face-on velocity dispersion 
of (Tv~150kms~^, indicating a highly turbulent gas environment. This is also suggested by 
the large R2i^35 ratio, expected in the high-pressure, highly turbulent ISM environments 
of strong mergers (Aalto ct al. 1995). Its CO SLED is highly excited up to CO J=4-3 
while the weak CO J=6-5 reported in Paper I is most likely due to pointing offsets as ZEUS 
measurements at the CSO indicate a much stronger line (Stacey 2011). 

No reasonable one-phase LVG fits were found for this source while attempts to fit the 
CO J=l-0 up to J=4-3 and R2i=35 yield highly unbound states (Kvii.~13-70). The best 
such fit: [Tkin,n,Kvir] = [70,10^13], yields r2i=0.96, 132=0.64, r43=0.32, the latter being -3x 
smaller than observed. The large of this system corresponds to Msf(H2)~6x 10^ Mq. On 
the other hand Xco~0.4Xi computed from the best LVG solution yields Mtot~5.2xl0^ Mq, 
and thus even the minimum mass expected for a dense and warm SF-fueling phase will 
dominate Mtot, as computed by a one- phase Xco factor. One- phase LVG solutions with their 
low densities, highly unbound dynamical states certainly do not hint this, while their poor 
fit of the global CO line ratios of this ULIRG suggest that at least two main gas phases are 
necessary to represent the average ISM state in this system. The presence of a massive gas 
phase much denser than those implied by one-phase LVG fits is further corroborated by an 
i"HCN/co=0.14 (HCN J=l-0 from Gracia-Carpio et al. 2008, CO 1-0 from Paper I) and an 
HCN(3-2)/(l-0) ratio of r32(HCN)=0.39 (Gracia-Carpio et al. 2008). 

Using only CO J=4-3, 3-2, and R2i=35 (the lower limit) to constrain the (h)-phase 
in a 2-phase model yields only warm Tkin=(90-150) K and dense n~(3x 10^-10^) cm^^ gas 
phases with Kvir~13-22, and Xco^~0. 85-1.40 (high values for the denser gas solutions). These 
solutions yield r^'^-^ j~2.3-3 for J-|-l=2,3,4 indicating CO lines with small/moderate optical 
depths, a result of the high temperatures and large Kvi,. values which, along with the high 
densities, are necessary for well-excited CO lines up to high-J levels and a high global R21 
ratio. The mass of such high-excitation gas can be computed from L^q^ =(l/(r43''))L43=(4- 
4.4)x109Li which yields Mh-cx(H2)~(3.4-6) x 10^ M©. Thus, as in some other (U)LIRGs 
studied here (e.g. IRAS 00057-1-4021), this galaxy contains large reservoirs of dense gas in 
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thermal and dynamical states much more extreme than those expected for dense gas in 
SF regions. As detailed in Paper I these could be indicative of other power sources such as 
turbulent heated regions (THRs) or CR-dominated regions (e.g. Papadopoulos 2010). Future 
high resolution high-J CO, ^^CO imaging of this system with ALMA would be very valuable 
for confirming and quantifying such mechanisms via determining the level of turbulence 
present in its disk (i.e. measuring crz(r) and the average CR energy density permeating its 
ISM. 

Despite the prominence of the dense highly-excited phase in the J=3-2, 4-3 lines, the 
(l)-phasc remains the main contributor of CO J=l-0 line luminosity in IRAS 17208-0014, 
with l['^1q~(8.7-9) X 10^ Li (~(66-68)% of the observed one). LVG fits of the corresponding 
"residual" ratio r2i~0.5, while Tkin-degenerate (as expected), all indicate low gas densities 
n~(l-3)xl0^cm~^ and a wide range of Kvir~l-70. The coresponding xS~(0.5-3)Xi (over 
the Kvir=l-20 range of LVG solutions) remains essentially unconstrained, with the lowest val- 
ues corresponding to the most dynamically unbound states as expected. Thus for Xco^=1.4 Xi 
(densest gas LVG solutions) we find Xco~'''^'*~(0.75-2.5) Xi. 

Using the states compatible with the HCN(3-2)/(l-0) ratio in IRAS 17208-0014 to de- 
fine the (h)-phase leads to even higher Xco~^''''~(2.5-5.9) Xi, with larger still Mh_ex (see Table 
1). These correspond to a high-density (~(3x 10^-10^) cm~^) phase with Kvir~l-6 (unphys- 
ical solutions with Kvir<l were not considered). Nevertheless, for Tkm=(75-105)K, LVG 
solutions with high-density gas (~10^cm^^) but Kvir=20 are also possible. These strongly 
unbound dense gas solutions yield xio-'~0.8Xi and similarly low Xco"^'^'' values (unless cold 
(l)-phase solutions corresponding to a SF-quiescent disk are considered). Observations of 
more HON lines and, crucially, of at least one H^^CN line can determine whether such 
extraordinary dynamical states are possible for the HCN- bright dense gas phase (see 3.4). 



B.25. IRAS 22491-1808 

This is a ULIRG with a spectacular morphology involving two nuclei separated by 
1.6" (~2.3kpc) (Carico et al. 1990), two tidal tails and numerous SF knots embedded in 
them (Scoville at el. 2000) indicating an advanced merger. It is also one of the very few 
ULIRGs where an AGN contributes more than half of the total IR luminosity (Farrah et al. 
2003). In our sample it stands out as the galaxy with the most highly excited CO J=3-2 
line with r32=1.88, though a large applied correction makes this value somewhat uncertain 
(see Paper 1). All one-phase LVG solutions correspond to dense gas n=(10*^-3xl0^) cm~^, 
with the best fits obtained for Tkin>40K. Unfortunately lack of more CO line ratios and 
especially at least one ^^CO line measurement permits a wide range of [Tkin, n, Kyir] states 



-58- 



up to Tkin=150 K and Kvir~l-20 to fit the observed ratios. The one-phase Xco values remain 
essentially unconstrained with Xco~(0.7-2.4) Xi and the low values associated with strongly 
unbound states (Kvir~13-22), while the low values with states that have Kvir~l-7. The 
corresponding molecular gas mass is Mtot~(6.3-22)xlO^M0 while Msf~5.2x10^ Mq. 

Another extraordinary aspect of IRAS 22491-1808 is that many of the conditions com- 
patible with its global CO line ratios indicate dense gas (~3xl0^) cm~^) with temperatures 
Tkm=(100-130) K, significantly larger than those of the dust (Tdust=(32-49) K for emissivity 
of /3=1.5, 2). Such potentially strong decoupling of gas and dust temperatures, especially 
at high gas densities where Tkin^Tjust is expected, are another indicator of other domi- 
nant power sources than PDR photons. Turbulent and CR heating are capable of driving 
large Tkin/Tdust>l inequalities even for the very dense gas mass reservoirs found in ULIRGs 
(Paper I, Papadopoulos 2010). These mechanisms can easily dominate when strong gas- 
rich/gas- rich spiral galaxy interactions drive the bulk of the molecular gas into compact 
regions in the two nuclei of the progenitors. This will inject large amounts of gas kinetic 
energy into ~100pc-sized regions fueling extreme turulence, while high SFR densities will 
establish extreme CRDRs. Interestingly two other double-nuclei systems, IRAS 12112-1-305 
and IRAS 08572-1-3915, also show very highly excited (i.e. rj_|_ij>l) CO lines. 

The IRAS 22491-1808 and similar systems are thus excellent cases for studying extreme 
ISM conditions and their drivers using ALMA. Such observations will be crucial for: a) infer- 
ring CR energy densities (by finding starburst region sizes), b) measuring velocity dispersions 
and the turbulence levels of molecular gas disks and c) directly providing key molecular line 
diagnostics for the presense of CR-dominated regions (CRDRs) and/or turbulent-heated 
regions (THRs) (Papadopoulos 2010; Bayet et al. 2011, Meijerink et al. 2011). 

B.26. NGC 7469 

This is a well-studied Sy 1 spiral galaxy whose AGN is surrounded by a nearly com- 
plete ring of powerful starburst activity contributing ~2/3 of its Ljr (Genzel et al. 1995; 
Riffel et al. 2006 and references therein). Interferometric CO J=l-0 imaging found its 
central starburst to be very gas- rich (Meixner et al. 1990; Tacconi & Genzel 1996). Its star- 
burst is thought triggered by an interaction with a neighboring galaxy IG 5283 that hes 83" 
(~27kpc) away (projected distance). A Xco~(l/5)Xco Qai was inferred for its inner ~lkpc, 
and attributed to the effects of strong SF feedback on molecular clouds (Genzel et al. 1995). 
A small Xco for the starburst region of the disk in NGC 7469 is indeed corroborated by 
CO J=3-2, 1-0 observations and 850/um/450/im dust continuum imaging (Papadopoulos & 
Allen 2000). However the same observations also revealed faint CO J=l-0 and 450/im dust 
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emission extending well beyond the staburst region, and containing a massive SF-quiescent, 
cold dust and gas reservoir with a Galactic Xco- 

These two ISM components make a good one-phase LVG fit impossible for the "cold" 
r2i=0.75 and the "warm" r43=0.83, r65=0.22 ratios observed. The best such solution with 
Tkin=70K, n(H2)=103cm-3, K,i,=A yields: 4^^^=0.97, 4^^^=0.75, ri3*^^=0.48, r^5^^=0.08, 
and rJ''^^=11. The first two are higher than those observed (~2o" higher in the case of T21), 
while for (4-3)/(l-0) and (6-5)/(l-0) the LVG-computed ratios are ~1.7 and ~2.75 times 
smaller than observed. The Xco=0.72 Xi value computed from the optimal one-phase LVG fit 
to the global CO SLED is dominated by the SF phase, and corresponds to Mtot~2.5x 10^ M0. 
The minimum SF gas mass expected on the other hand is Msf=6.5x 10^ M©. 

Using only the CO J=3-2, 4-3, G-si and the ^^CO J=2-l lines to constrain the 
high-excitation (h)-phase in a 2-phase models yields a typical LVG solution for Tkin=60 K, 
n(H2)=10'^ cm~^, and Kvir=13 with a corresponding Xco'*=0.66 Xi. From this solution we 
compute l[5'(^'i_o=2.51x109Li, and L^J^ i_o=0.97x 10^ Li which yield p^!r^^~0.39. Then us- 
ing the Galactic Xco=5Xi deduced for the extended CO J=l-0 and submm continuum dust 
emission of the disk in NGC7469 (Papadopoulos & Allen 2000), we obtain xio"^^^~1.88 Xi, 
which is ~2.6x higher than that obtained by the one-phase LVG fit of the global SLED. 

Thus the cold extended disk of NGC7469, despite containing ~3x more molecular gas 
mass than its central starburst, is nearly incospicious in the global CO SLED, resulting to a 
significant underestimate of the total molecular gas mass in this galaxy. Only the available 
spatial information about its presence could rectify this, a state of affairs mirrored in other 
such systems (e.g. NGC1068). For disk-dominated LIRGs with strong excitation gradients 
induced by the presense of central starbursts only CO, ^'^CO hne and dust continuum imaging 
of their ISM can avoid such pitfals (see 3.4). This is quite unlike ULIRGs where it is the low- 
excitation "cloud-envelope" -phase that contains the smaller mass fraction and is concomitant 
or closely follows the distribution of a massive much denser gas component. 



B.27. IRAS 23365+3604 

Near-IR imaging of this LIRG reveals tidal tails indicating a merger, and a single nucleus 
embedded into a nearly face-on disk that is nearly ~20kpc (17.3") in diameter, and shows 
disturbed spiral structure (Surace et al. 2000). The nucleus is classified as a LINER (Veilleux 
et al. 1995). The gas dynamics as revealed by CO J=l-0 interferometry (DS98) is strongly 



■^The wide beam of the CSO was used to obtain the CO 6-5 measurent for this extended object 
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decoupled from stellar dynamics (Genzel et al. 2001), another hallmark of strong mergers. 
Its molecular gas disk is compact with ~1.2kpc diameter, and a face-on velocity dispersion 
of -lOOkms-i (DS98). Its CO SLED has a "cool" CO(2-1)/(1-0) ratio (=0.75) but a well- 
excited CO J=3-2 line with a "warm" CO(3-2)/(l-0) (=0.82), while its large R2i(=17) is 
typical for mergers. One-phase LVG models yield moderately good but rather degenerate 
fits with [Tkin,n,Kvir] = [30-65, 10^, 13] and [70-150,3x10^, 2.2] being the optimal solution 
ranges. The best solutions are obtained for Tkin=50K, 70K with corresponding line ratios: 
4^^)^0.88, 0.84, rf2''^^=0.56,0.55, ri3^^=0.24, 0.27 and ^^=20, 16. These are compatible 
with the measured values though the models typically yield "warmer" r2i and "cooler" r32 
ratios than observed, in the latter case by ~l.lcr. The corresponding Xco~(l-1.3) Xi, which 
gives Mtot~(7.3-9.5)xlO'^Mo, while Msf~2.5x109 M©. 

The latter represents a substantial fraction of the total molecular gas mass where much 
higher gas densities must be prevailing. A clear indication of another highly-excited denser 
gas phase is provided by CO J=6-5 line which is ~9-16 times more luminous than expected 
from the best LVG solutions. The presence of a dense gas phase able to emit such high-J 
CO lines is also indicated by a high HCN/CO J=1^0 line ratio of -I^jjcn/co^ 

0.10 (HCN from 

Gracio-Carpio et al. 2008, CO from Paper I). Using the CO J=3-2, J=4-3 (upper limit), 
J=6-5 and ^^CO J=2-l lines to constrain the high-excitation phase yields: [Tkin, n, Kvir] = [40- 
80, 10^, 40] as the optimum solution range, with Tkin=60K being the best. The corresponding 
CO SLED has r^i^=1.94, rj^ = 1.90, rJ3^=1.50, and r[;5^=0.45, typical of low/modest optical 
depths (indicatively Tio~0.2, T32~2.5), a result of the exceptionally large K^ir. A less optimal 
solution range with lower densities and more modest (but still high) Kyir values does exists 
[Tkin, n, Kvir] = [90-150, 3x10^, 22], but could not reproduce the high HCN(l-0)/CO(6-5) ratio 
of Rhcn/co65~0.48±0.18 observed for this galaxy as such a low-density phase would hardly 
emit any HCN J=l-0 (Rhcn/co65<0.08). 

The mass of the dense, warm, and strongly unbound phase can be found from X^o^~0.45 Xi 
(for Tkin=60 K), and l[5^'i_o=(1/42^)Lco,3-2~3.15x 10^ Lj, which yields Mh_cx~l-4xl09 Mq. 
This large gas mass with such high K^ij. and high densities is rather surprising since such 
strongly unbound states are associated with cloud-envelopes or intercloud diffuse molecular 
gas whose n<10^cm~^. The remaining CO J=l-0 luminosity of the (l)-phase on the other 
hand is Lqq ;^_q~4. 15x10^ Li with an essentially unconstrained Xco. Setting x[iQ=(0.5-5) Xi 
to encompass the range from a warm, diffuse and unbound gas component to an underlying 
cold, self-gravitating, disk phase yields X^o"^^^~(0.5-3) Xi and Mtot~(3.7-22) xlO'^ M©. The 
low (and thus Mtot) values are more probable because of the negligible CO J=2-l lumi- 
nosity of the (l)-phase (a cold, self-gravitating, phase would have higher J=2-l luminosities). 

Using the SNR-GMC interfaces as benchmark entities for dense yet highly unbound 
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(due to the kinetic energy injected by SNR shocks) molecular gas phase makes clear that a 
mass fraction of Mh_cx/Mtot~0.06-0.38 exceeds that expected per GMC for a SNR-impacted 
gas phase of a CMC's mass). It is worth noting that only in the powerful radio galaxy 

3C293 galaxy-sized molecular gas reservoirs at high densities were inferred to be at such 
strongly gravitationally unbound states, likely caused by the powerful jet-lSM interaction. 
While no such ACN influence on the ISM is apparent in IRAS 23365+3604 it is interesting to 
note that it is one of the 3 ULIRGs (the others being Arp 220 and Mrk 273) where molecular 
gas velocity fields decoupled from those of stars and ionized gas are found, with molecular 
gas having ~2 times larger linewidths (Cohna et al. 2005). High- J CO, ^^CO and HCN and 
H^^CN line observations are necessary for confirming such extraordinary dynamical states 
for the massive dense gas reservoir residing in IRAS 23365-1-3604 and similar systems. 
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Fig. 1. — The Xco factors, estimated from Equations 3 and 5, using the results from one- 
phase LVG radiative transfer models. The black-bar histogram corresponds to the best and 
least degenerate solutions obtained from these models. 
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Fig. 2. — The distribution of the iyn(H2)/Tr_i_o factor for the LVG solutions obtained for 
the sample (black bars as in Figure 1). 
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Fig. 3. — The Xco factors, computed from Equation 3 and the 1-phase LVG solutions obtained 
for the CO hne ratios measured in LIRGs (see Paper I) for two different dynamical regimes 
of average gas motions: self-gravitating (1 < Kvir<4), and unbound (Kvir>4). 
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Fig. 4. — The Xco factor, estimated from Equation 6 using the /3 j+i j values computed from 
one-phase radiative transfer models. 
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Fig. 5. — The minimum molecular gas masses necessary for fueling and Eddington-limited 
star formation in the LIRGs of our sample (see section 2.3). 




Fig. 6. — The gas mass fraction expected at overdensities x>Xo where x=n/ (n) (see Equation 
13) in supersonic turbulent gas, from spiral disks (M=10) up to ULIRGs (M=300). 



Table 1. Molecular gas masses and Xco factors 



Name 


-"-^^tOt V CO J 


Mprp''^(x^o"P''^)'^ 




IV'Tt., 


R emarks 




(10^ Mn) 


(10^ Mn) 


(10^ Mn) 


(10^ Mn) 




IRAS 00057+4021 


4.5 (1.1) 


3.4-8.2(0.83-2.0) 


1.2 


1.6 


degenerate xS values 


IZwl 


8.6 (1.5) 


9.5 (1.65) 


0.4 


2.8 


Galactic xS adopted 


NGC 828 


4.6-6.3 (0.8-1.1) 


4-14 (0.7-2.5) 


0.3 


0.3 


degenerate Xco 


IRAS 02483+4302 


1 6 (0 45) 


4.3-8.6 (1.2-2.4) 


1.4 


3.3 


degenerate X^l^ values 


IRAS 03359+1523 


22 (2.5) 




0.8 




SF-auiescent ISM 


VII Zw 031 


14.5 (1.25) 


39 (3.4) 


2.0 


3.1 


Cold extended disk? 


IRAS 05189-2524 


13.7 (3.5) 


9.4 (2.4) 


3.6 


3.6 


(li)-pliase SLED assumed i 


IRAS 08030+5243 


14.3-42.9 (1.5-4.5) 




2.4 




degenerate 1-phasc solutions ^ 


IRAS 08572+3915 


4.8-9.6 (3-6) 




3.9 




Highly excited CO J=6-5 line 


Arp 55 


2.9-9.2 (0.25-0.8) 




0.9 




degenerate l-phase solutions 


UGC 05101 


1.9-2.9 (0.4-0.6) 




1.6 






NGC 3310 


0.03-0.15 (0.4-2.2) 




0.1 




degenerate Xpo values 


IRAS 10565+2448 


3.8-4.8 (0.6-0.75) 


4.8 (0.75) 


2.6 


2.6 


(h)-pliase SLED assumed 


Arp 299 


1.2 (0.42) 




1.9 




Unconstrained Xco 


IRAS 12112+0305 


15.2-53.5 (1.5-5.3) 


7.1-34.3 (0.7-3.4) 


4.8 


4 


Mrk231 


1.8-3.2 (0.25-0.45) 


21-49 (3-7) 


5.5 


8-25 


CO 6-5, HCN lines used in 2-phase model 


Arp 193 


1.6-3.5 (0.35-0.75) 


3.3-8.1 (0.7-1.8) 


0.9 


1.7-7.1 


CO 3-2, HCN lines used in 2-phase model 


NGC 5135 


1.1-1.4 (0.35-0.45) 




0.3 






Mrk273 


10.4 (2.0) 


13 (2.5) 


3.3 


6 


Cold extended disk? 


3C293 


2.3-7.4 (0.5-1.6) 


13.4-18 (2.9-3.9) 


0.05 


1.1-3.5 


Galactic xS adopted 


IRAS 14348-1447 


11-17 (0.65-1) 


24-30 (1.4-1.7) 


5.6 


5.6 


degenerate (l)-phase 


Zw 049.057 


0.30-0.52 (0.35-0.6) 




0.47 







Table 1 — Continued 



Name 




TV /r /"v \ a 
Mtot(Xco) 

(iU Mq) 


T,^(2-ph)/Y(2-ph)vb 
Mtot (Xco ') 

/"I n9 i\/r ^ 


TV /r c 
MsF 

/'in9i\/r ^ 


TV /r d 
Mh-ex 

/"I n9 i\/r ^ 


Remarks 


Arp 220 




1.85 (0.30) 


14.7-27.5 (2.4-4.5) 


4 


12.6-25.2 


HON lines used in 2-phase model 


NGC6240 




2.5 (0.30) 


8.4-27.8 (1-3.3) 


1.5 


4.9-24.5 


HCN lines used in 2-phase model 


IRAS 17208- 


-0014 


5.2 (0.40) 


9.8-33 (0.75-2.5) 


6.2 


3.5-6 


degenerate Xco values 


u u 






34-77 (2.6-5.9) 


u 


27-70 


HCN lines used in 2-phase modeP 


IRAS 22491- 


-1808 


6.3-22 (0.70-2.4) 




5.2 




degenerate Xco values 


NGC7469 




2.5 (0.72) 


6.5 (1.88) 


0.65 


1.65 


Cold extended disk present 


IRAS 23365+3604 


7.3-9.5 (1-1.3) 


3.7-22 (0.5-3) 


2.5 


1.4 


Unconstrained xS 



^Total molecular gas mass from Equation 3 and the best one-phase LVG model parameters (the corresponding Xco 
value in Xi units), with average values adopted in cases of significant LVG solution range degeneracy. 

''Total molecular gas mass from Equation 10 and a 2-phase fit (the corresponding Xco value in Xi units). 

■^The minimum molecular gas mass necessary for an Eddington-limited star formation rate (see 2.3). 

'^The gas mass of the (h)-phase in a 2-phase model when this is used. If a standard (h)-phase SLED and mass 
normalization is used (see 2.4) then Mh-cx=MsF- 

*^For LVG solutions with Kvir(HCN)-l-6 (see B.24). 



